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Most  natural  fluid  bodies,  the  ocean,  the  atmoshpere,  many  lakes 
and  most  estuaries,  exhibit  inherent  stratification.  Recently,  with 
the  rapid  development  of  coastal  and  offshore  resources  and  with  the 
recognition  of  environmental  quality  as  a  fundamental  consideration, 
stratified  flows  have  taken  on  a  new  importance  for  engineers  and 
hydrologists. 

In  this  dissertation,  the  mixing  region  between  two  parallel  and 
homogeneous  streams  of  fluid  having  different  densities  and  velocities 
is  investigated  experimentally.  The  flow  is  produced  in  a  flume  with  a 
test  section  24.38  m  long,  1.22  m  high  and  0.61  m  wide  by  pumping  the 
upper  0.61  m  deep  layer  in  a  closed  circuit.  The  density  of  each  layer 
is  controlled  by  the  addition  of  sodium  chloride  in  such  a  manner  that 
the  lower  layer  is  about  2.5  percent  denser  than  the  upper  layer. 


XI 


Tests  are  conducted  at  two  downstream  locations  and  at  several 
velocities  between  2.7  cm/sec  and  7.8  cm/sec.  The  mean  velocity 
and  the  fluctuating  velocity  components  are  measured  with  hot 
film  anemometers  while  the  density  structure  is  measured  with  a 
high  resolution  conductivity  probe.  Profiles  of  mean  velocity, 
viscous  stress,  and  Reynolds  stress  are  obtained  through  the  mixing 
region.  Energy  spectra  of  the  fluctuating  velocity  components 
are  also  obtained  at  several  points  within  the  region  and  under 
several  flow  conditions. 

Over  the  range  of  conditions  investigated,  the  interface  moves 
about  and  distorts  a  small  amount  in  a  temporally  and  spatially 
random  fashion.  However,  it  is  always  distinct  and  readily 
identifiable;  there  are  no  large  convoluting  motions  of  the  inter- 
face. The  velocity  structure  is  similar  in  some  ways  to  flow 
over  a  smooth  rigid  surface. 

The  results  show  that  for  the  conditions  investigated  the 
interfacial  boundary  layer  is  marginally  stable.  When  a  well- 
developed  turbulent  boundary  layer  exists  above  the  interface, 
the  Richardson  number  at  the  inflection  point  in  the  velocity 
profile  is  of  order  one  or  less. 
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CHAPTER  I 
INTRODUCTION 

Motivation 

Studies  of  stratified  flows  are,  with  few  exceptions,  motivated 
by  two  important  facts.  First,  most  natural  fluid  bodies  are,  to  some 
extent,  stratified.  The  ocean,  the  atmosphere,  many  lakes,  and  most 
estuaries  exhibit  inherent  stratification;  ostensibly,  stratification 
is  the  rule,  not  the  exception.  Secondly,  when  effluents  are  dis- 
charged into  the  water  or  air  there  is  generally  a  difference  in  den- 
sity between  the  effluent  and  the  receiving  fluid;  this  produces  a 
stratified  flow. 

Most  aspects  of  oceanic  and  atmospheric  circulation  are  intimate- 
ly tied  to  the  density  structure  of  the  fluid.  For  this  reason,  the 
dynamics  of  stratified  flows  have,  for  many  years,  been  important  to 
meteorologists  and  oceanographers.  More  recently,  with  the  rapid 


Stratification  denotes  a  layered  structure.  There  may  be  a  finite 
number  of  distinct  layers  in  which  case  the  system  is  referred  to  as 
an  "N-layer"  system.  In  the  limit,  the  fluid  may  have  an  infinite 
number  of  infinitely  thin  layers  in  which  case  the  term  "continuously 
stratified"  might  be  used.  The  property  differentiating  between  lay- 
ers can  in  general  be  any  property  of  the  fluid.  However,  in  this 
paper,  "stratified"  will  imply  density  stratification,  since  the  den- 
sity is  the  ultimate  factor  affecting  the  dynamics  of  the  flow.  The 
density  variations  can,  in  turn,  result  from  temperature  variations 
and/or  variations  in  dissolved  or  suspended  constituents.  Nothing 
has  been  said  regarding  the  orientation  or  extent  of  the  layers  with 
respect  to  other  layers  or  the  gravitational  field;  in  this  paper, 
the  term  is  used  loosely  to  imply  a  nonhomogeneous  density  structure. 


development  of  coastal  and  offshore  resources  and  with  the  recognition 
of  environmental  quality  as  a  fundamental  consideration,  the  subject 
has  taken  on  a  new  importance  for  engineers  and  hydrologists. 

Many  of  the  world's  large  population  centers  are  located  on  coast- 
lines near  estuaries  or  on  lakes.  These  centers  depend  on  their 
neighboring  hydraulic  systems  for  transportation,  food,  waste  disposal, 
recreation,  and  a  multitude  of  other  societal  requisites.  Many  times 
these  requirements  are  mutually  exclusive  when  the  system  is  improperly 
managed.  The  cornerstone  of  proper  management  is  an  accurate  and  reli- 
able model  of  the  flow  within  the  lake,  estuary,  or  coastal  zone. 
Nevertheless,  many  mathematical  models  used  today  are  unable  to  accu- 
rately describe  the  flows  when  stratification  is  present.  In  some 
models,  this  is  due  to  the  use  of  vertically  averaged  equations  of 
motion.  These  equations  cannot  account  for  changes  in  flow  direction 
or  speed  within  the  fluid  column.  However,  even  those  models  that 
attempt  to  account  for  stratification  by  employing  a  multi-layer  ap- 
proach suffer  from  an  incomplete  knowledge  concerning  the  transport  of 
mass,  momentum,  and  energy  between  layers  in  a  stratified  fluid.  Conse- 
quently, the  mathematical  models  available  today  require  extensive 
calibration  and  may  be  unreliable  if  significant  variations  are  made 
from  the  calibration  condition.  If  the  dynamics  of  natural  fluid  sys- 
tems are  to  be  understood,  modeled  effectively  and  reliably,  the  effect 
of  stratification  must  be  understood. 

The  consequences  of  stratification  are  generally  three-fold: 
1.  The  mean  flow  pattern  is  altered  by  the  redistribution  of 
inertia  and  pressure. 


2.  Turbulent  energy  is  enhanced  or  decreased,  depending  on  the 
nature  of  the  stratification.  This,  in  turn,  alters  the 
transfer  of  mass,  momentum,  and  energy. 

3.  Stratification  allows  for  the  existence  of  gravity  waves 
internal  to  the  fluid.  These  internal  waves  are  capable  of 
transmitting  energy  and  momentum  across  large  distances  with 
little  mixing  or  dissipation  of  energy. 

Any  one  or  all  three  of  the  above  can  have  a  profound  effect  on  the 
dynamics  of  the  fluid  system.  For  an  extensive  discussion  of  the  many 
phenomena  associated  with  stratified  fluids,  the  reader  is  referred  to 
Turner  (1973),  Phillips  (1969),  and  Yih  (1965).  These  books  all  have 
lengthy  bibliographies  for  further  reference. 

With  the  increasing  impact  man  exerts  upon  nature  and  therefore 
himself  comes  the  need  to  better  understand  the  dynamics  of  natural 
fluid  systems  and  the  effects  of  density  stratification. 

Two-Layer  Flows 

Many  studies  of  stratified  flows,  including  this  one,  are  con- 
cerned with  various  aspects  of  two-layer  flows.  This  type  of  flow, 
which  is  characterized  by  the  presence  of  two  thick  layers  of  nearly 
homogeneous  fluid  separated  by  a  thin  transition  region  in  which  the 
density  changes  abruptly,  continues  to  receive  considerable  attention 
for  a  number  of  important  reasons. 

The  first  and  most  basic  reason  is  the  geometric  simplicity  of 
the  flow.  As  with  all  investigations  of  complicated  phenomena,  strat- 
ified flows  must  first  be  broken  down  into  basic  structural  units; 


when  these  units  are  understood,  more  complicated  phenomena  can  be 
considered. 

Secondly,  and  perhaps  more  important  in  terms  of  the  engineering 
application,  the  two-layer  flow  closely  represents  many  natural  and 
man-made  occurrences.  For  example,  salt-wedge  and  fjord-type  estuaries 
often  exhibit  a  density  structure  which  is  closely  modeled  as  a  two- 
layer  flow.  Additionally,  when  heated  v/ater  or  effluents  from  indus- 
trial processes  are  discharged  into  the  environment  the  resulting  flow 
many  times  approximates  a  two-layer  flow,  the  dynamics  being  control- 
led in  large  part  by  the  processes  taking  place  at  the  interface. 
Finally,  the  formation  and  evolutionary  development  of  the  oceanic 
thermocline  and  the  atmospheric  inversion  are  strongly  dependent  on 
the  mixing  process  at  the  thermocline  or  inversion.  These  phenomena 
are  often  modeled  as  two-layer  stratified  flows.  This  is,  of  course, 
a  limited  list  of  the  possible  occurrences;  many  other  examples  exist. 

A  third  reason  for  studying  two-layer  flows,  which  has  far- 
reaching  implications  for  the  study  of  continuously  stratified  fluids, 
is  the  observance  of  step-like  density  structures  in  many  continuously 
stratified  flows.  Using  high  resolution  instruments  to  document  the 
fine  scale  structure  of  continuously  stratified  fluids,  several  inves- 
tigators have  found  that  in  many  natural  fluid  bodies  the  vertical 
density  and  velocity  profiles  are  not  smoothly  varying,  as  was  origi- 
nally expected,  but  are  step-like  in  appearance.  Thick  layers  of 
nearly  homogeneous  fluid  are  separated  by  thin  transition  regions  in 
which  the  velocity  and  density  change  abruptly.   In  the  ocean  this 


The  term  "fine  structure"  is  sometimes  used  in  the  literature  when 
referring  to  this  type  of  stratification. 


structure  has  been  documented  by  Woods  (1968),  Woods  and  Wiley  (1972), 
and  Osborne  and  Cox  (1972).  Woods  referred  to  the  nearly  homogeneous 
regions  as  "layers"  and  the  high  gradient  regions  as  "sheets."^  This 
fine  structure  has  also  been  observed  numerous  times  in  the  atmosphere; 
see,  for  example.  Browning  (1971),  Browning  and  Watkins  (1970),  and 
Emmanuel  (1973).  The  step-like  structure  has  even  been  observed  in 

the  fluff  layer  of  a  tidal  channel  where  stratification  is  a  conse- 

2 

quence  of  fine  suspended  sediments.   It  is  increasingly  apparent  that 

this  density  structure  is  a  widespread  phenomenon,  and  the  formation 
depending  little,  if  at  all,  on  the  stratifying  agent.  This  step- 
like density  structure  might  possibly  be  modeled  as  the  superposition 
of  many  two-layer  flows. 

A  Regime  Concept  for  Two-Layer  Flows 

In  fluid  mechanics,  the  concept  of  regimes  is  a  useful  tool  for 
describing  or  defining  the  characteristics  of  a  flow.  For  example, 
flows  are  often  described  as  falling  within  the  compressible  or  in- 
compressible, laminar  or  turbulent  regimes.  In  the  study  of  two-layer 
stratified  shear  flows,  a  differentiation  can  be  made  on  the  basis  of 
interfacial  dynamics.  Observations  reported  by  numerous  investigators 
Thorpe  (1968),  Koop  (1976),  Kantha  (1975),  and  Browand  and  Wang  (1972), 
to  mention  just  a  few,  reveal  that  interfacial  instabilities  can  pro- 
duce two  distinctly  different  modes  of  interfacial  mixing. 


In  this  paper  the  high  gradient  region  may  be  alternatively  referred 
to  as  the  interfacial  region. 

2 
Seminar  presented  at  the  University  of  Florida  by  Dr.  W.  R.  Parker 

of  the  Institute  of  Oceanographic  Sciences,  November  18,  1977. 


First,  under  conditions  of  high  shear  and/or  low  static  stability, 
large  vortical  structures,  known  as  Kelvin-Helmholtz  billows,  form  at 
the  interface.  These  large  coherent  structures  produce  an  overturning 
of  the  density  interface,  entraining  fluid  from  the  upper  and  lower 
layers  into  the  vortex  spiral.  Many  studies,  beginning  with  Helmholtz 
(1868)  and  Kelvin  (1871)  and  continuing  to  the  present,  have  dealt 
with  various  aspects  of  this  instability.  Recent  experimental  studies 
in  both  liquids  and  gases--Thorpe  (1973,  1969),  Roshko  (1976),  Koop 
(1976),  Browand  and  Winant  (1973),  Browand  and  Wang  (1972),  Scotti  and 
Corcos  (1969)--have  greatly  increased  the  understanding  of  the  initia- 
tion, growth,  and  decay  of  this  type  of  instability.  Some  of  these 
papers  have  excellent  photographic  plates  showing  the  Kelvin-Helmholtz 
billow.  The  correlation  between  experiment  and  theory  strongly  sub- 
stantiates the  idea  that  these  structures  are  the  result  of  inviscid 
instabilities.  In  this  regime,  the  interface,  in  an  unbounded  flow, 
diffuses  rapidly  with  time  and/or  space  until  at  some  point  it  again 
becomes  stable. 

The  second  regime  exists  when  the  shear  is  small  and/or  the  static 
stability  is  large.  In  this  case  the  mixing  seems  to  result  from  the 
breaking  of  small  wavelets  which  develop  a  sharp  cusp-like  shape  at  the 
crest  and/or  trough,  ejecting  fluid  from  one  layer  into  the  other.  The 
interface  in  this  regime,  unlike  the  first  regime,  remains  distinct  and 
easily  identifiable.  The  mixing  process,  in  this  case,  more  closely 
approximates  turbulent  entrainment  than  diffusion.  The  origin  and 
nature  of  the  motions  which  produce  this  mode  of  mixing  are  not  as  well 
understood  as  those  which  produce  the  Kelvin-Helmholtz  billow.  They 
might  be  the  end  result  of  viscous  or  inviscid  instabilites  which  can 


originate  at  the  interface  or  at  a  distant  location.  Benjamin  (1963) 
has  identified  three  types  of  instabilities  which  can  arise  in  the 
region  of  a  flexible  boundary.  His  results  indicate  that  viscous 
instabilities  of  the  Tollmein-Schlichting  type  can  arise  before  the 
Kelvin-Helmholtz  instability.  Browand  and  Wang  (1973),  however,  have 
proposed  that  the  two  different  types  of  mixing  result  from  two  dif- 
ferent modes  of  an  inviscid  instability. 

Thus,  because  of  the  dissimilarities  between  the  two  modes  of 
mixing  and  the  fundamental  questions  which  still  remain,  it  is  useful 
at  this  time  to  study  and  model  the  two  regimes  separately.  Eventual- 
ly, as  more  data  are  obtained,  it  may  be  possible  to  consider  the  two 
regimes  through  one  comprehensive  model . 

Scope  of  This  Investigation 

The  intent  of  this  dissertation  is  to  document  experimentally,  as 
accurately  as  possible,  the  density  and  velocity  structure  of  the  inter- 
facial  region  in  a  weakly  turbulent  stratified  shear  flow.  Previous 
investigations,  many  of  which  are  discussed  in  Chapter  II,  have  either 
provided  qualitative  descriptions  of  the  mixing  process,  documented 
interlayer  transfer  rates  for  mass  momentum  and  energy  or  established 
stability  parameters  for  transition  from  laminar  to  turbulent  flow. 
Few  have  accurately  measured  the  density  and  velocity  structure  in  the 
transition  region  between  the  layers.  This  study  documents  the  time 
averaged  mean  density  and  velocity  profiles,  the  turbulent  and  viscous 
shear  stress  profiles  and  the  turbulent  energy  spectra  through  the 


In  this  paper,  when  the  terms  "regime  one"  or  "regime  two"  are  used 
alone,  the  reference  will  be  to  the  conditions  described  above. 
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interfacial  region  of  a  weakly  turbulent  (regime  2)  stratified  shear 
flow.  The  intent  is  to  provide  accurate  baseline  data  needed  to 
establish  sound  theories  describing  the  entrainment  of  fluid  from  one 
layer  to  another  or  the  shear  stress  at  the  interface  between  fluids 
with  different  densities  and  velocities.  These  two  generic  problems 
form  the  basis  for  more  specific  problems  such  as  the  dynamics  of 
salt-wedge  and  fjord-type  estuaries,  the  development  and  evolution  of 
the  oceanic  thermocline  and  the  dynamics  of  atmospheric  inversions. 

In  Chapter  II,  many  papers  published  to  date  are  categorized  and 
several  important  results  from  these  studies  are  briefly  discussed. 
Chapter  III  describes  the  facilities,  instrumentation,  and  procedures 
used  to  obtain  the  data  presented  in  Chapter  IV.  The  last  chapter 
discusses  the  conclusions  reached  from  this  study  and  makes  recommenda- 
tions regarding  future  work  in  this  area. 


CHAPTER  II 
BACKGROUND  DISCUSSION 

Introduction 

Numerous  articles  on  the  dynamics  of  flow  near  a  density  interface 
have  been  published  in  a  wide  variety  of  journals.  Generally,  these 
studies  fall  into  one  of  three  categories:  (1)  studies  of  laminar 
flow,  (2)  investigations  into  the  dynamic  stability  of  two-layer  strat- 
ified shear  flows,  and  (3)  studies  concerned  with  the  turbulent  trans- 
fer of  mass,  momentum,  and  energy  across  a  density  interface.  This 
last  group  can  be  further  divided  on  the  basis  of  where  the  turbulent 
energy  originates.  Sometimes  mixing  is  induced  by  a  turbulent 
energy  source  located  far  from  the  density  interface.  This  situation 
is  usually  characterized  by  the  absence  of  mean  flow  vorticity  in  the 
region  of  the  interface.  Other  times  the  transfer  results  from  the 
production  of  turbulent  energy  at  the  interface.  In  this  case,  a 
mean  velocity  gradient  is  required  at  the  interface  in  order  to  main- 
tain production  of  turbulence.  Still  other  times  both  energy  sources 
are  present.  This  is  the  most  common  condition  found  in  nature; 
however,  it  is  also  the  most  complicated. 

Studies  of  laminar  flow  and  hydrodynamic  stability  are  based  on 
firmly  established  and  highly  developed  mathematical  theories.  The 
theory  of  turbulent  flow,  however,  has  not  yet  been  developed  to  the 
point  where  complex  flows  can  be  predicted  from  basic  equations. 
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Investigations  of  turbulent  flow,  therefore,  must  rely  on  experimental 
observation  augmented  by  less  well  developed  mathematical  models. 

Although  there  is  no  direct  link  between  linear  stability  theory 
and  fully  developed  turbulent  flow,  much  insight  into  the  origin  of 
turbulent  motions  can  be  obtained  by  examining  the  results  of  linear 
stability  theory.  This  is  the  primary  reason  for  including,  in  this 
chapter,  the  brief  discussion  on  stability. 

Laminar  Flow 

Keulegan  (1944)  and  Lock  (1951)  analyzed  the  laminar  boundary 
layer  between  two  streams  of  fluid  with  differing  density  and  viscos- 
ity. The  fluids  were  assumed  to  be  incompressible,  nondiffusive  and 
inmiscible,  such  that  the  density  and  viscosity  were  discontinuous 
across  the  interface.  They  used  the  standard  boundary  layer  assump- 
tions and  solved  for  the  velocity  in  each  layer  by  matching  the  solu- 
tions at  the  interface.  As  one  would  expect,  their  solutions  show 
that  the  boundary  layer  grows  with  the  square  root  of  the  downstream 
distance,  the  constant  of  proportionality  depending  on  the  viscosity 
in  the  respective  layer.  The  exact  shape  of  the  velocity  profile  is 
dependent  on  the  product  of  the  density  and  viscosity  ratios.  Since 
the  requirement  of  laminar  flow  is  seldom  realized,  studies  of  this 
type  have  limited  direct  practical  application  toward  problems  in 
engineering  and  geophysics.  They  are  important,  however,  since  they 
supply  the  basic  mean  flow  required  for  stability  analysis,  and  they 
form  a  limiting  case  upon  which  to  build  more  applicable  models. 
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stability  of  Stratified  Shear  Flows 

Helmholtz  (1868)  and  Kelvin  (1871),  interested  in  the  generation 
of  water  waves  by  wind,  were  the  first  to  investigate  analytically  the 
stability  of  stratified  shear  flows.  Their  studies  examined  the  case 
where  the  density  and  velocity  are  discontinuous  at  the  interface  and 
uniform  within  each  layer  (the  stratified  vortex  sheet).  Later,  Taylor 
(1931),  Goldstein  (1931),  and  Holmboe  (1962),  to  list  some  of  the  more 
widely  quoted  examples,  analyzed  piecewise  linear  velocity  and  density 
profiles.  Holmboe  (1962)  and  Miles  (1961,  1963)  extended  the  analysis 
to  include  hyperbolic  tangent  velocity  and  density  profiles.  For  a 
lengthy  discussion  on  the  stability  of  parallel  inviscid  flows  the 
reader  is  referred  to  the  review  article  by  Drazin  and  Howard  (1966). 

Two  important  general  observations  grow  out  of  these  analytic 
solutions.  The  first  is  Miles'  1/4  Theorem:  if  the  gradient  Richardson 
number  is  everywhere  greater  than  1/4,  the  flow  will  be  stable.  This 
is  only  a  sufficient  condition  for  stability;  it  states  nothing  about 
when  the  flow  will  become  unstable  or  if  stability  will  persist  below 
1/4.  The  second  observation  is  more  qualitative  in  nature.  The  point 
at  which  the  flow  becomes  unstable  and  the  nature  of  the  instability 
depend  strongly  on  the  details  of  the  velocity  and  density  structure 
or,  alternatively,  the  Richardson  number  profile.  This  second  point  is 
reinforced  by  the  results  from  experimental  studies  of  two-layer 
stratified  flows.  Figure  1  shows  two  hypothetical  velocity  and  density 
profiles  and  their  respective  Richardson  number  profiles.  Although  the 
total  change  in  velocity  and  density  is  the  same  for  both  cases,  the 
shapes  of  the  Richardson  number  profiles  are  radically  different.  This 
difference  results  solely  from  the  diversities  of  the  vertical  length 
scales. 
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Thorpe  (1969,  1973),  Scotti  and  Corcos  (1969),  and  Delisi  and 
Corcos  (1973)  investigated  experimentally  the  stability  of  flows  where 
the  gradient  Richardson  number  is  a  minimum  at  the  interface  (see 
Figure  1).  They  found  that  Miles  Theorem  held  and  that  Kelvin- 
Helmholtz  billows  formed  when  the  interface  became  unstable.  The  in- 
vestigations of  Browand  and  Wang  (1972)  and  Wang  (1975),  however, 
reveal  a  much  different  stability  boundary  in  the  Richardson  number-- 
wave  number  plane.  It  should  be  noted  that  in  these  experiments  the 
viscous  boundary  layer  was  much  thicker  than  the  diffusion  layer. 
Under  these  conditions,  the  profile  of  the  gradient  Richardson  number 
will  exhibit  a  local,  or  possibly  an  absolute,  maximum  at  the  inter- 
face (see  Figure  1).  Koop  (1976),  using  the  same  facility  as  Wang 
(1975),  demonstrated  that  two  distinctly  different  modes  of  instability 
form  at  the  interface,  depending  on  the  initial  conditions.  These  two 
instabilities  correlate  well  with  the  results  of  Hazel's  numerical 
study  (1972)  and  the  analytic  study  of  Holmboe  (1962).  The  character 
of  the  interfacial  disturbance  created  by  these  two  instabilities 
corresponds,  at  least  qualitatively,  to  the  two  regimes  discussed  in 
Chapter  I. 

To  this  point,  only  inviscid  stability  theory  has  been  discussed; 
the  effects  of  viscosity  need  to  be  considered.  Maslowe  and  Thompson 
(1971)  and  Gage  (1973)  numerically  analyzed  the  stability  of  stratified 
shear  flows,  including  the  effects  of  both  viscosity  and  heat  conduc- 
tion. Their  results  support  the  idea  that  viscosity  has  a  stabilizing 
influence.  It  should  be  noted,  however,  that  the  velocity  and  density 
profiles  that  were  used  again  led  to  a  Richardson  number  profile  with 
a  minimum  at  the  interface. 
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Benjamin  (1963)  investigated  the  stability  of  a  flexible  boundary 
in  an  inviscid  flow,  and  extended  his  results  to  the  case  of  an 
inviscid  fluid  flowing  over  a  viscous  fluid.  His  model  assumed  a 
uniform  density  in  both  the  upper  and  lower  layers,  with  a  Blasius 
velocity  profile  in  the  upper  layer  and  zero  flow  below  the  interface. 
The  Richardson  number  profile  under  these  assumptions  is  undefined  at 
and  below  the  interface.  This  model,  nevertheless,  closely  approxi- 
mates many  actual  laboratory  and  geophysical  flows.  When  the  density 
difference  is  large,  as  with  air  over  water,  Benjamin's  results  show 
that  the  Kelvin-Helmholtz  (K-H)  instability  will  be  the  dominant  in- 
stability. However,  when  the  density  difference  is  small,  Tollmien- 
Schlichting  (T-S)  instabilities  can  arise  at  half  the  velocity  required 
for  a  K-H  instability.  Since  the  T-S  instability  requires  energy  dis- 
sipation to  grow,  viscosity  must  be  considered  a  destabilizing  factor 
under  certain  conditions.  To  date,  no  one  to  the  author's  knowledge 
has  examined  analytically  or  numerically  the  stability  of  an  unbounded 
viscous  stratified  shear  flow  where  the  Richardson  number  profile  has  a 
maximum  at  the  interface.  It  may  be  that  under  this  condition,  viscos- 
ity is  a  destabilizing  factor,  and  T-S  instabilities  can  arise  before 
K-H  instabilities.  This  could  explain  the  two  regimes  discussed  in 
Chapter  I:  regime  one  is  undoubtedly  the  result  of  K-H  instabilities; 
regime  two  might  be  the  result  of  T-S  instabilities.  As  Benjamin  (1963) 
demonstrates,  the  K-H  instability  extracts  energy  from  the  mean  flow  at 
a  much  higher  rate  than  the  T-S  instability;  therefore,  if  the  K-H  in- 
stability exists,  it  will  dominate. 

Keulegan  (1949)  experimentally  investigated  the  stability  of  a 
two-layer  viscous  stratified  flow.  Using  an  approach  which  paralleled 
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Jeffreys'  theory  (1926)  for  the  generation  of  water  waves,  Keulegan 

1 


derived  a  stability  parameter, 

(vg  Ap/p) 


1/3 


6 


U 

where  v  is  the  viscosity,  g  the  gravitational  acceleration,  p  the  fluid 

density,  Ap  the  density  difference  across  the  interface,  and  U  the 

velocity  of  the  upper  layer.  His  experimental  data  showed  that  the 

critical  value  for  the  generation  of  interfacial  wavelets  is: 

e  =  0.127,  Re  <  450 
e  =  0.178,  Re  >  450 

where  Re  is  the  Reynolds  number  based  on  the  hydraulic  radius  of  the 

flume.  The  instabilities  observed  by  Keulegan  had  regime  two 

characteristics. 

Turbulent  Entrainment  Across  Density  Interfaces 

Many  investigators,  beginning  with  Rouse  and  Dodu  (1955)  and 
Cromwell  (1960),  have  studied,  using  oscillating  grid  tanks,  the 
problem  of  turbulent  mass  transport  across  density  interfaces.  Within 
these  tanks,  turbulent  motions  are  generated  above  and/or  below  a 
stable  density  interface  by  oscillating  a  grid  of  rigid  bars.  This 
experimental  apparatus  makes  the  measurement  of  entrainment  rates  rel- 
atively simple;  however,  since  no  mean  velocity  is  present,  measurement 
of  turbulent  properties  is  difficult.  Hopfinger  and  Toly  (1976)  were 
able,  nevertheless,  to  measure  the  turbulent  fluctuations  using  a  hot- 
film  anemometer  mounted  on  a  rotating  arm.  Their  results  show  that  the 


This  parameter  is  sometimes  referred  to  as  the  Keulegan  number  (K) 
where  K  =  e^. 
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turbulent  structure  is  closely  tied  to  the  amplitude  and  frequency  of 
oscillation  and  to  the  geometry  of  the  grid.  The  turbulent  energy  is 
concentrated  in  frequency  bands  close  to  the  grid  oscillation  frequency 
and  harmonics  of  this  frequency.  Dispersion  of  the  turbulent  energy  to 
wider  frequency  bands  is  slow  in  this  type  of  facility  since  there  is 
no  mean  velocity  shear,  and  therefore  only  limited  vortex  stretching. 
Turbulent  motions  produced  in  this  manner  can  generate  standing  waves 
at  the  interface  which  grow  and  break  due  to  partial  resonance  with 
the  oscillator. 

Despite  these  drawbacks  some  significant  findings  have  resulted 
from  these  experiments.  Linden  (1975)  found  that  the  rate  of  increase 
in  the  potential  energy  of  the  fluid  column  is  proportional  to  the  rate 
at  which  turbulent  energy  is  supplied  to  the  interfacial  region.  There 
is  no  simple  relationship,  however,  between  the  rate  of  change  of  po- 
tential energy  and  the  energy  input  at  the  grid.  He  also  found  that  as 
much  as  half  the  energy  supplied  to  the  interface  could  be  radiated 
away  by  internal  waves  if  the  Vaisala  frequency  below  the  interface  is 
not  zero. 

Turner  (1968)  and  Crapper  and  Linden  (1974)  studied  the  effect  of 
molecular  diffusivity  on  the  mixing  rate.  They  found,  for  low  turbu- 
lent Reynolds  numbers, 

Re^  E  ;^  <  50  , 

V 

that  the  mechanics  of  the  mixing  process  and  the  rate  of  mixing  were 
dependent  on  the  turbulent  Peclet  number: 


Pe 


Ui£ 
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where  U, ,  is  a  turbulent  velocity  scale,  Ji  is  a  turbulent  length 
scale,  k  the  coefficient  of  molecular  diffusivity,  and  v  the  kinematic 
viscosity.  When  Pe  <  200,  the  interface  has  a  diffussive  core  and 
the  thickness  of  the  interface  depends  on  the  value  of  Peclet  number. 
When  Pe  >  200,  there  is  no  dependence  on  Peclet  number;  the 
turbulent  eddies  completely  penetrate  the  interface,  and  molecular 
diffusion  across  the  interface  is  insignificant.  Turner  (1968)  also 
found  when  Peclet  number  is  large  (the  case  when  salt  is  used  as  a 
stratifying  agent)  that: 

'  "A 

where  Ue  is  the  entrainment  velocity  and  the  Richardson  number  is 
based  on  U, ,  and  l. 

Ellison  and  Turner  (1959),  Kato  and  Phillips  (1969),  and  Moore 
and  Long  (1971)  have  investigated  the  problems  of  turbulent  entrainment 
across  a  density  interface  in  the  presence  of  a  mean  velocity  shear. 
Their  results  support  the  relationship: 

Ue   p.-l 

The  discrepancy  between  the  above  two  power  laws  has  been  the  source 
of  much  debate  [see  Long  (1975)  ].  It  is  not  clear  which  expression 
is  correct  or  if,  as  Long  (1975)  argues,  the  difference  is  due  to  the 
presence  of  a  mean  velocity  shear.  It  is  clear,  however,  that,  in 
general , 

j^  =   f(Ri,  Pe,  Re)  and 
"l 
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that  only  when  Reynolds  number  and  Peclet  numbers  are  large  is 

ir=f(Ri)  . 

Dynamic  Structure  of  a  Stratified  Mixing  Layer 

To  date,  there  is  no  model  for  the  turbulent  stratified  half-jet 
or  mixing  layer  that  is  valid  over  a  wide  range  of  Richardson  numbers. 
There  are,  however,  two  limiting  cases  that  have  been  studied  and 
modeled  extensively.  First,  as  the  Richardson  number  approaches  zero, 
the  homogenous  turbulent  half- jet  results.  At  the  other  extreme,  the 
Richardson  number  approaches  infinity  and  the  static  stability  of  the 
interface  becomes  so  large  that  the  interface  acts  much  like  a  solid 
plate.  The  flow  of  air  over  water  is  an  example  of  this  extreme  case. 
The  turbulent  structures  exhibited  by  these  two  limiting  flows  are 
vastly  different.  For  example,  the  eddy  viscosity  is  almost  constant 
near  the  center  of  the  homogeneous  half-jet,  while  at  a  rigid  plate 
the  eddy  viscosity  approaches  the  fluid  viscosity.  When  a  turbulent 
stratified  two-layer  flow  is  in  regime  one,  its  dynamic  structure  is 
more  closely  approximated  by  the  homogeneous  half-jet  than  flow  over 
a  flat  plate;  however,  when  in  regime  two  the  stiuation  is  reversed. 
This  dissertation  is  concerned  with  the  dynamic  structure  of  a  regime 
two  stratified  mixing  layer.  Therefore,  the  many  aspects  of  a  regime 
one  mixing  layer  will  not  be  discussed.  Readers  interested  in  this 
type  of  flow  are  referred  to  Roshko  (1976)  and  Brown  and  Roshko  (1974) 

Lofquist  (1960)  made  extensive  measurements  of  velocity  and  den- 
sity through  the  interfacial  region  of  a  regime  two  stratified  flow. 
Because  hot  film  anemometers  were  in  their  infancy  at  the  time  his 
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work  was  done,  Lofquist  was  able  to  directly  measure  only  the  pressure 
gradient  and  the  mean  velocity  and  density  profiles.  The  interfacial 
shear  stress  was  calculated  using  the  equations  of  motion  in  conjunc- 
tion with  the  above  measurements.  All  of  Lofquist's  experiments 
reveal  a  velocity  structure  that  closely  resembles  the  flow  over  a 
flat  plate. 

Csanady  (1978)  demonstrated,  using  Lofquist's  (1960)  data,  that 
a  "Law  of  the  Wall"  could  be  formulated  for  density  interfaces.  The 
"Law  of  the  Interface,"  however,  differs  from  the  traditional  "Law  of 
the  Wall"  in  that  it  depends  on  two  parameters,  Reynolds  number  and 
Keulegan  number.  Csanady  chose  to  formulate  the  "Law  of  the  Inter- 
face" as, 

where  U^  is  the  interfacial  friction  velocity,  <  the  Von  Karman  con- 
stant, z  the  vertical  coordinate,  and  v  an  effective  viscosity 
defined  by, 

2 
V, 


'e   dU 
dz 


interface 


B(K)  and  v  (K)  are  constants  which  depend  on  the  turbulent  Keulegan 
number  defined  by, 

Apgv 
K  = 


pU. 


3 


where  v  is  the  kinematic  viscosity,  Ap  the  density  difference  across 
the  interface,  and  p  the  fluid  density. 
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Dimensional  Analysis 

In  this  section,  a  functional  form  of  the  velocity  profile  in  the 
region  of  a  density  interface  is  developed  from  dimensional  reasoning. 
The  approach  parallels  the  classical  analysis  of  turbulent  flow  over  a 
rigid  surface,  except  that  the  roughness  length  used  in  the  analysis 
of  a  rigid  surface  is  replaced  by  a  buoyancy  length  scale. 

Under  the  Boussinesq  approximation,  and  assuming  the  velocity 
profile  in  the  region  of  the  density  interface  is  only  a  function  of 
the  fluid  properties,  the  shear  stress  at  the  interface,  gravity  and 
the  vertical  coordinate,  the  following  holds: 

U  =  f(Pj^,  y,  g',  U^,  d 
where : 

p  is  the  fluid  density  of  the  lower  layer; 

y  is  the  fluid  viscosity  (assumed  to  be  equal  in  upper  and 

lower  layers); 

g'  is  the  reduced  gravity  =  gAp/p^; 

g  is  the  acceleration  of  gravity; 

Ap  is  the  density  jump  across  the  interface; 

U*  is  the  interfacial  friction  velocity  =  "^^f^i^ 

T.  is  the  interfacial  shear  stress; 
1 

5  is  the  vertical  distance  from  the  interface;  and 
U  is  the  mean  velocity. 

The  above  independent  variables  (excluding  c)  can  be  combined  to 
form  two  length  scales:  a  viscous  length  scale. 


'V  "  P„U 
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and  a  buoyancy  length  scale. 


A  physical  interpretation  of  the  buoyancy  length  scale  is  as 
follows:  z.    is  the  vertical  distance  an  inviscid  particle  of  density 
p  would  travel,  in  a  medium  of  density  p^-Ap,  if  it  was  given  an 
initial  vertical  velocity  U^.  Therefore,  if  U^  is  a  measure  of  the 
turbulent  velocity  scales,  £.  is  a  measure  of  the  vertical  distortions 
of  the  interface  caused  by  the  turbulent  motions  within  the  flow. 

Paralleling  the  approach  used  in  the  analysis  of  flow  over  a 
rigid  surface,  two  alternative  functional  expressions  for  the  mean 
velocity  profile  near  a  density  interface  are: 
U(0 


and 


U* 


=  f 


where 


and 


^-"l"  is  a  Richardson  number, 
is  a  Reynolds  number. 


I     I 

^■1  3  =  -r~  ~  ^   ""S  ^   Keulegan  number. 

^r*   ^b 


When  K  »  1,  the  viscous  length  scale  dominates  the  buoyancy  length 
scale.  Therefore,  the  distortions  of  the  density  interface  (the 
interfacial  wavelets)  will  be  deeply  imbedded  within  a  viscous  sublayer 
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and  should  have  little  influence  on  the  flow.  Under  this  condition, 
the  first  functional  expression  would  appear  most  appropriate.  When 
K  <<  1 ,  the  distortions  of  the  interface  would  be  much  larger  than  the 
viscous  sublayer  and  would  therefore  dominate  the  flow.  Under  this 
condition,  the  second  expression  would  appear  most  appropriate. 

The  results  from  Lofquist  (1960)  and  from  this  study  (Figure  30) 
indicate  the  existence  of  a  region  above  the  interface  where  the  veloc- 
ity profile  has  a  logarithmic  shape.  It  should  therefore  be  possible 
to  represent  the  velocity  profile  in  this  region  by: 

Ulil=  A(K)  £n  (  '-^   )  +  B(K)  ,  K  >  1  . 

In  the  case  of  flow  over  a  rigid  surface, 

A(K)  =  1/k 
where  k  is  Von  Karman's  constant,  there  is  no  reason  to  assume,  in  the 
case  of  flow  over  a  density  interface,  that  A  is  constant.  The  flux 
of  mass  through  the  interface  could  alter  the  turbulent  structure  and 
therefore  the  value  of  A. 

The  above  expression  is  similar  in  form  to  the  "Law  of  the  Inter- 
face" given  by  Csanady  (1978),  except  that  Csanady  assumed  A  =  l/< 
and  replaced  the  fluid  viscosity  with  an  "equivalent  viscosity."  The 
author  feels  the  above  expression  is  a  more  physically  meaningful  form. 
The  functions  A(K)  and  B(K),  as  well  as  the  range  of  e.   over  which  this 
expression  holds,  remain  to  be  determined..  Finally,  it  seems  reason- 
able to  assume  that  as  K  becomes  large,  A(K)«=^1/k. 

Governing  Equations 

The  governing  equations  for  turbulent  stratified  flow  are  derived 
from  the  differential  form  of  the  continuity  and  Navier-Stokes  equations 
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using  Reynolds  decomposition  technique.  Specific  concern  here  is  for 
two-layer  flow  in  a  retangular  channel  where  the  mean  velocity  is  nearly 
unidirectional.  The  x-axis  is  positive  in  the  primary  flow  direction 
while  the  y-axis  is  positive  upward  and  the  z-axis  is  perpendicular  to 
the  X  and  y  axes  using  a  righthanded  coordinate  system. 

If  the  flow  is  assumed  to  be  nearly  two  dimensional  and  free  of 
secondary  currents  then  at  the  center  of  the  tank  (z=0),  where  the  mea- 
surements presented  in  Chapter  IV  were  made,  the  continuity  equation  is 

^  +  ^=0  2.1 

3x   ay 

where  U  and  V  are  the  time  averaged  mean  velocities  in  the  x  and  y 
directions,  respectively. 

If  the  pressure  throughout  the  fluid  is  assumed  to  be  hydrostatic, 
the  y-component  momentum  equation  is  simply  the  hydrostatic  equation 

^  =   -g  Po(x,y),  2.2 

where  Po  and  po  are  the  time  mean  pressure  and  density,  respectively,  and 
g  is  the  acceleration  of  gravity. 

Finally,  if  third  order  correlations  and  correlations  between  veloc- 
ity and  density  fluctuations  are  assumed  small,  then  the  mean  advective 
accelerations,  the  horizontal  pressure  gradient,  the  side  wall  shear 
stress  and  the  interfacial  shear  stress  are  the  only  significant  terms 
in  the  x-component  momentum  equation.  Furthermore,  if  po  is  approxi- 
mately equal  to  p^,  the  Boussinesq  approximation  holds  and  the  x-compon- 
ent momentum  equation  is 

9X    3y    p^  3x    3y     3Z   '  ^-^ 
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where 

is  the  shear  stress  maintained  by  the  side  walls  and 

xy         Pu  ^^  ' 
is  the  shear  stress  maintained  by  the  interface. 

Equation  2.3  can  be  further  simplified  if  one  assumes  that  the  flow 
does  not  vary  with  distance  downstream  and  the  side  wall  stresses  are 
insignificant.  Under  these  conditions  equation  2.3  becomes 

Pu  9x    3y   • 

Appendix  E  contains  the  full  form  of  the  governing  equations  and 
presents  in  detail  the  derivation  of  equations  2.1  through  2.4.  In  addi- 
tion, a  seven  parameter  expression  for  the  shear  stress  profile  at  the 
center  of  the  tank  (z=0)  is  developed  from  equations  2.1  through  2.3 
assuming  a  self-similar  velocity  and  density  profile.  This  expression 
for  the  shear  stress  will  be  discussed  further  in  Chapter  IV. 
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FIGURE  1:  THE  EFFECT  OF  VERTICAL  LENGTH  SCALES  ON 
THE  RICHARDSON  NUMBER  PROFILE 


CHAPTER  III 
EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Facilities 

The  internal  wave  tank  at  the  University  of  Florida's  Coastal  and 
Oceanographic  Engineering  Laboratory  was  used  in  this  study.  The  tank 
was  designed  as  a  multi-purpose  stratified  flow  facility;  consequently, 
it  had  available  many  operational  features  which  were  not  needed  for 
this  experiment.  Only  those  aspects  of  the  facility  pertinent  to  the 
execution  of  this  study  will  be  discussed  here.  For  a  detailed  de- 
scription of  all  the  capabilities  of  this  tank,  see  Sheppard,  Shemdin, 
and  Wang  (1973). 

Figure  2,  a  schematic  drawing  of  the  overall  facility,  shows  the 
locations  of  the  main  components.  The  test  section  is  24.38  meters 
long,  1.83  meters  high,  and  0.61  meter  wide  in  the  lower  two-thirds  of 
the  tank,  while  the  upper  one-third  is  0.91  meter  in  width.  Figure  3 
shows  a  cross-sectional  view  of  the  test  section.  In  this  study,  the 
bottom  one-third  is  filled  with  salt  water,  and  the  middle  one-third 
with  fresh  water.  The  upper  one-third  of  the  tank  is  a  wind  tunnel 
designed  for  the  generation  of  wind  waves  and  wind-induced  currents; 
it  was  not  used  during  this  experiment.  The  sides  of  the  test  section 
are  constructed  from  one-half  inch  (1.27  cm)  thick  glass  panels,  allow- 
ing visual  observation  of  phenomena  from  both  sides  of  the  tank. 
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The  flow  is  produced  by  pumping  the  upper  layer  in  a  closed  cir- 
cuit. The  pumping  system  consists  of  two  centrifugal  pumps  connected 
in  parallel.  The  60  hp  pump  can  produce,  in  the  test  section,  veloci- 
ties approximating  0.5  m/s.  However,  at  low  flow  rates  (velocities 
less  than  8  cm/s)  this  pump  heats  the  water  excessively  and  is  diffi- 
cult to  control.  The  second  pump  is  only  5  hp  but  it  can  produce  flows 
up  to  8  cm/s;  therefore,  for  low  flow  rates  the  small  pump  was  used  to 
minimize  thermal  stratification  in  the  tank. 

The  saline  solution  is  mixed  before  each  experiment  in  one  of 
two  4,000  gallon  (15,140  liter)  storage  tanks;  the  fresh  water  is 
stored  in  the  other.  The  saline  solution  is  usually  dyed  with  red 
food  color,  so  that  the  interface  between  the  fresh  and  salt  water  is 
easily  identified.  From  the  storage  tanks  the  water  is  fed  by  gravity 
into  the  internal  wave  tank.  The  salt  water  enters  the  tank  through  a 
horizontal  slot  in  the  bottom;  its  inflow  rate  can  be  regulated  from 
0  to  5  gal/min  (0.32  liter/sec)  to  an  accuracy  of  +  0.05  gal/min 
(+  .0032  liter/sec)  using  a  Fischer  and  Porter  tapered  tube  flow 
meter.  Fresh  water  from  the  city  water  supply  can  be  continuously 
added  to  the  upper  layer  through  a  tap  in  the  12  inch  (30.48  cm)  pipe 
which  is  located  upstream  of  the  pumps.  The  free  surface  level  is 
maintained  by  the  discharge  of  fresh  v/ater  across  a  weir  at  the  down- 
stream end  of  the  test  section.  This  continuous  exchange  of  fresh 
water  forestalls  any  significant  change  in  the  upper  layer  density  due 
to  the  entrainment  of  salt  water  across  the  density  interface.  The 
maximum  fresh  water  inflow  rate  is  limited  to  12  gallons/min 
(0.76  liter/sec). 
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The  entrance  to  the  test  section  has  been  significantly  modified 
for  this  experiment  from  that  described  by  Sheppard  et  a1 .  (1973). 
The  new  entrance  section  was  designed  with  two  criteria  in  mind: 
(1)  to  allow  for  a  smooth  transition  in  the  flow  where  the  fresh  upper 
layer  meets  the  lower  salt  layer,  and  (2)  to  produce  a  velocity  pro- 
file which  is  straight,  parallel,  and  as  close  as  possible  to  the 
fully  developed  profile  for  flow  in  a  rectangular  channel.  The  intent 
of  the  second  criterion  is  to  minimize  longitudinal  variations  in  the 
velocity  profile.  Figure  4  shows  the  general  layout  of  the  entrance 
section.  The  profile  shaping  element  is  constructed  from  880  pieces 
of  3/4-inch  (1.91-cm)  I.D.  PVC  pipe.  Each  tube  is  cut  to  a  prescribed 
length  based  on  its  position  within  the  tank  cross-section.  For  de- 
tails of  the  construction  see  Appendix  A, 

The  transition  plate,  which  is  just  downstream  of  the  profile 
shaping  element,  is  horizontal  at  their  point  of  contact;  the  down- 
ward slope,  as  well  as  the  curvature,  increases  with  distance  down- 
stream. This  design  was  employed  to  insure  that  the  upper  and  lower 
layers  meet  smoothly  and  with  no  significant  mixing,  even  though  the 
vertical  position  of  the  interface  may  vary  several  centimeters  due  to 
set-up  and  entrainment.  The  end  of  the  transition  plate  is  always  at 
least  25  cm  below  the  interface. 

The  exit  from  the  test  section  has  also  been  modified  for  this 
study.  Figure  5  shows  the  layout  at  the  downstream  end  of  the  tank. 
This  design  minimizes  saltwater  entrainment  and  blockage  close  to  the 
interface  while  separating  the  upper  and  lower  layers.  By  rotating 
the  spliter  wedge  about  its  pivotal  point,  the  vertex  of  the  wedge  can 
be  corredpondingly  adjusted  to  match  long-term  changes  in  the  vertical 
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position  of  the  interface.  The  vertically  slotted  weir,  located  10  to 
20  cm  downstream  of  the  wedge  vertex,  is  designed  to  enhance  turbulent 
mixing.  This  eliminates  the  build-up  of  a  wedge  of  intermediate  den- 
sity water  at  the  downstream  end  of  the  tank,  thus  reducing  blockage 
effects  near  the  interface. 

Instrumentation 

The  instrument  system  for  this  investigation  consists  of  six  in- 
dependent transducers,  nine  signal  conditioners,  a  traverse  mechanism, 
and  an  analog-to-digital  (A/D)  tape  recorder.  The  output  signal  from 
each  transducer  is  passed  through  one  or  more  signal  conditions  and  is 
sampled  and  recorded  onto  magnetic  tape  by  the  Kinemetrics  Model  DDS- 
1103  A/D  tape  recorder.  Figure  6  shows  the  general  arrangement  of  the 
instrument  package. 

The  Kinemetrics  A/D  recorder  can  simultaneously  sample  up  to 
16  channels  of  data  at  any  one  of  75  predetermined  sample  rates  between 
,05  and  1,000  samples/sec.  The  input  voltage,  which  must  be  between 
-5.0v  and  +5.0v,  is  resolved  into  4,096  discrete  values.  Consequently, 
the  input  resolution  for  this  unit  is  limited  to  +0,0012  volts. 

The  traverse  mechanism  is  driven  by  a  variable  speed  electric 
motor  which  allows  the  speed  to  be  continuously  controlled  between 
0.02  cm/sec  and  2.0  cm/sec.  It  can  move,  vertically,  a  total  distance 
of  71.0  cm.  The  probes  for  transducers  two  through  six,  which  measure 
temperature,  conductivity,  and  velocity,  are  mounted  at  the  end  of  the 
traverse  support  rod. 

The  first  transducer,  which  indicates  the  vertical  position  of 
the  traverse  support  rod,  consists  simply  of  a  10-turn,  lOk-ohm 
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variable  resistor  in  series  with  two  fixed  resistors;  the  potentiameter 
is  attached  to  the  gear  drive  of  the  traverse  mechanism.  The  unit  has 
a  linearity  of  0.05  percent  over  the  range  of  the  output  voltage  from 
-5.0v  to  +5.0v.  The  signal  from  the  position  indicator  goes  to  a  low- 
pass  filter  with  20-Hz  cutoff  frequency.  The  output  from  the  filter 
is  then  input  to  the  Kinemetrics  tape  recorder.  Since  this  voltage 
can  be  determined  only  to  within  +0.0012  volt  by  the  A/D  recorder,  the 
resolution  of  the  position  indicator  is  limited  to  +0.085  mm. 

The  second  transducer  is  a  thermistor;  its  electronic  schematic 
diagram  is  shown  in  Figure  7.  The  probe.  Model  SF2-42A68  manufactured 
by  Victory  Engineering  Company,  has  a  sensor  only  1.5  mm  in  diameter 
and  a  time  constant  of  0.4  second  in  still  water.  The  thermistor  was 
calibrated  against  a  mercury  thermometer  with  a  resolution  of 
+0.05  degrees  Fahrenheit.  The  thermistor  resolution  is  assumed,  there- 
fore, to  be  the  same.  A  typical  calibration  curve  for  this  transducer 
is  shown  in  Figure  8,  and  although  the  electronic  circuit  was  designed 
to  produce  a  linear  output,  a  small  degree  of  nonlinearity  can  be  seen 
in  the  calibration  curve.  To  compensate  for  this,  the  calibration 
points  were  fitted  with  a  second  order  polynomial  (three  coefficients), 
using  the  method  of  least  square  error.  The  difference  between  any 
one  of  the  calibration  points  and  the  curve  fit  value  was  always  less 
than  the  resolution  stated  above.  The  output  from  the  thermistor 
passes  to  a  low-pass  filter  with  a  cutoff  frequency  of  20  Hz;  the  out- 
put from  the  filter  is  recorded  by  the  A/D  recorder. 

The  electrical  conductivity  of  the  fluid  is  measured  by  a  third 
transducer.  Knowing  the  conductivity  and  temperature,  the  density  of 
the  fluid  can  be  determined.  The  conductivity  probe,  shown  in 
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Figure  9,  is  constructed  of  acrylic  plastic  and  has  one  stainless  steel 
and  one  platinum  electrode.  When  the  probe  is  operating,  the  surround- 
ing fluid  is  continuously  drawn  into  the  interior  chamber  through  a 
0.5  mm  slit  and  a  0.3  mm  diameter  vertical  hole  at  the  bottom  of  the 
probe.  The  fluid  exits  through  the  plastic  tube  connecting  the  probe 
and  the  vacuum  pump.  The  electrical  current  density  between  the  elec- 
trodes is  small  except  in  the  vertical  hole,  where  it  becomes  wery 
large.  Thus,  the  resistance  between  the  electrodes  at  any  instant  is 
governed  by  the  conductivity  of  the  fluid  in  the  vertical  hole,  not  by 
the  resistance  at  the  fluid-electrode  interface.  This  reduces  the 
effects  caused  by  electrode  contamination  and  plating.  Figure  10  shows 
the  electronic  circuit  diagram  for  the  transducer.  The  horizontal  slit 
in  the  probe  is  designed  to  enhance  the  selected  withdrawl*  of  the 
surrounding  fluid.  For  more  on  the  characteristics  of  this  unit  see 
Sheppard  and  Doddington  (1977). 

The  major  advantage  of  this  design  over  the  so-called  "single 
electrode"  conductivity  probe  used  by  previous  investigators  is  that 
the  D.C.  drift  problem  has  been  virtually  eliminated  while  still  main- 
taining the  high  vertical  resolution  needed  for  the  study  of  density 
stratified  fluids.  There  are,  however,  two  important  disadvantages. 
First,  the  time  response  for  this  probe  is  much  slower  than  the  single- 
point  conductivity  probe  because  the  fluid  must  be  physically  drawn 
through  the  hole  in  the  probe.  Based  on  a  comparison  of  temperature 
and  conductivity  time  series  data,  a  drop-off  in  the  response  curve  is 


♦When  a  sink  is  placed  in  a  stably  stratified  fluid,  it  tends  to  draw 
fluid  only  from  its  own  level  provided  the  velocity  is  not  too  large. 
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believed  to  occur  at  about  2  Hz.  Second,  the  probe  is  physically  larg- 
er than  a  single-point  conductivity  probe;  therefore,  its  presence  will 
alter  the  flow-field  for  a  greater  distance.  In  turn,  this  necessi- 
tates a  greater  spacing  between  the  hot  film  probes  and  the  conductiv- 
ity probe.  As  a  consequence  of  both  the  above,  correlations  between 
density  and  velocity  fluctuations  measured  with  this  system  are  of 
limited  significance. 

A  calibration  curve  of  conductivity  voltage  versus  density  for  a 
fluid  temperature  of  75.0  degrees  Fahrenheit  is  shown  in  Figure  11, 
while  Figure  12  shows  the  effect  of  temperature  on  conductivity  voltage 

for  two  different  levels  of  salt  concentration.  Using  these  two  plots, 

_5 
the  density  can  be  determined  to  within  +5.0  x  10   gm/cc. 

The  output  voltage  from  this  transducer  is  again  input  to  a  low- 
pass  filter  with  a  20-Hz  cutoff  frequency,  and  the  output  from  the 
filter  is  again  recorded  on  magnetic  tape. 

The  last  three  transducers  consist  of  three  hot  film  sensors  and 
their  respective  anemometers.  The  first  two  sensors  are  part  of  an 
x-type  hot  film  probe,  Thermo-Systems  Model  1243-60w  NaCl .  This  probe 
is  used  to  measure  both  the  mean  velocity  and  the  fluctuating  compon- 
ents of  velocity  in  the  plane  of  the  sensors.  The  third  hot  film. 
Model  1232w  NaCl,  is  a  wedge-type  probe  used  only  for  mean  velocity 
measurements.  All  three  sensors  use  Thermo-Systems  Model  1050  general 
purpose  anemometers. 

The  output  from  each  anemometer  is  first  passed  to  a  Thermo- 
Systems  Model  1057  signal  conditioner.  This  conditioner  suppresses 
the  excessive  D.C.  Voltage  to  make  the  signal  compatible  with  the  in- 
put requirements  of  the  A/D  recorder.  From  the  conditioner,  the 
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signal  is  sent  to  a  low-pass  filter  with  a  cutoff  frequency  of  17  Hz 
and  from  there  the  signal  goes  to  the  A/D  recorder  where  it  is  digi- 
tized and  recorded  onto  magnetic  tape. 

The  purpose  of  the  low-pass  filters  is  to  prevent  contamination 
of  the  signal  by  high  frequency  noise,  specifically  60-cycle  noise, 
which  is  picked  up  from  the  large  pumps,  electric  motors,  and  other 
support  equipment  in  the  laboratory.  The  low-pass  filters  used  with 
the  position,  temperature,  and  conductivity  transducers  are  fourth- 
order  Butterworth  filters  with  the  -3  db  point  at  20.2  Hz.  At  10  Hz, 
the  power  loss  through  these  filters  is  less  than  two  percent.  Sixth- 
order  Butterworth  filters  are  used  with  the  hot  film  anemometers.  The 
-3  db  point  for  these  filters  is  at  17.3  Hz,  which  assures  that  the 
power  loss  at  10  Hz  is  less  than  one  percent.  Spectral  analysis  of 
the  velocity  data  indicates  that  for  these  experiments,  negligible 
energy  exists  above  7  Hz;  therefore,  these  low-pass  filters  have  little 
or  no  effect  on  the  true  signal  from  their  respective  transducers. 

Data  Acquisition  Procedure 

Two  different  techniques  were  employed  to  obtain  data  during  these 
experiments.  The  first  is  referred  to  as  a  constant  speed  traverse 
(CST),  while  the  second  is  called  a  stepped  traverse  (ST). 

The  primary  objective  of  the  constant  speed  traverse  was  to 
obtain,  in  a  quasi-instantaneous  manner,  the  structure  of  the  density 
transition  region.  This  was  accomplished  by  moving  the  probe  package 
through  the  interface  at  a  nearly  constant  speed  between  0.5  mm/sec 
and  1.5  mm/sec.  As  the  probes  moved  down,  the  position,  conductivity, 
and  thermistor  voltages  were  each  sampled  and  recorded  at  the  rate  of 
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25  Hz.  Since  the  thickness  of  the  interfacial  region  was  about  1.5  cm, 
it  took  no  more  than  30  seconds  to  traverse  the  entire  transition  re- 
gion. The  second  objective  of  the  CST  was  to  determine  the  position 
of  the  interface,  defined  as  the  point  where  the  density  equals  the 
average  density  between  the  upper  and  lower  layers.  The  position  of 
the  interface  was  then  plotted  as  a  function  of  time  and  all  other 
measurements  were  referenced  to  the  interface. 

The  stepped  traverse  was  used  to  obtain  information  on  the  struc- 
ture of  the  turbulent  boundary  layer  in  the  region  of  the  density 
interface.  With  this  method  the  probes  were  placed  at  a  pre-determined 
position  and  the  time  series  data  from  the  thermistor,  the  conductivity 
probe,  and  the  three  hot  film  anemometers  were  recorded.  The  time 
series  were  sampled  at  the  rate  of  15  Hz  over  a  period  of  540  seconds-- 
a  total  of  8,100  samples.  While  these  data  were  being  recorded  on 
magnetic  tape,  the  position  voltage  was  recorded  in  the  lab  book  for  a 
later  determination  of  the  exact  probe  position.  The  probes  were  then 
moved  to  another  position  and  the  time  series  data  were  again  recorded. 
This  procedure  was  repeated  until  sufficient  points  were  obtained 
(15  to  20)  to  determine  the  profile  through  the  transition  region. 

The  drop  in  the  position  of  the  interface  over  the  540  seconds 
required  to  record  the  time  series  data  at  each  point  was  much  less 
than  0.2  mm.  This  is  insignificant  when  determining  the  position  of 
the  probe  relative  to  the  interface.  However,  the  total  time  required 
to  complete  an  ST  profile  was  as  long  as  three  hours;  care  was  taken 
to  account  for  the  drop  in  the  interface  over  this  time  interval. 
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The  steps  followed  in  taking  data  were: 

1.  The  salt  water  was  allowed  to  flow  in  under  the  fresh  water 
at  the  rate  of  3.5  gal/min;  at  this  rate  15  hours  were  re- 
quired to  fill  the  tank.  The  filling  process  began  the  night 
before  the  tests. 

2.  The  instruments  were  calibrated  and  mounted  on  the  traverse 
just  prior  to  the  start  of  the  experiment. 

3.  The  pump  was  started  and  the  flow  increased  to  the  velocity 
desired.  This  required  as  much  as  one  hour  since  it  had  to 
be  done  slowly  to  prevent  surges  in  the  tank  and  seiching  of 
the  interface. 

4.  The  pump  was  run  for  an  additional  hour  or  more  before  any 
data  were  taken.  This  allowed  time  for  any  transients  to  die 
out. 

5.  A  CST  was  taken  first,  followed  by  an  ST  profile,  which  was 
followed  in  turn  by  another  CST.  This  sequence  of  three 
profiles  was  usually  repeated  several  times;  each  ST  was 
preceded  and  followed  by  a  CST. 

Data  Reduction  Procedure 

The  data  sets  were  reduced  on  the  IBM  System/370  computer  at  the 
University  of  Florida  using  the  Fortran  IV  language.  The  first  step  in 
the  reduction  process  was  to  convert  the  instrument  voltages  into  phys- 
ically significant  values.  Because  of  nonlinearities  in  most  of  the 
calibration  curves  (the  position  indicator  is  the  only  instrument  which 
is  linear  over  the  entire  range  of  operation),  this  must  be  done  before 
any  other  numerical  computations.  To  make  this  conversion,  the  least 
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squares  method  was  used  to  fit  an  N   order  polynomial  to  the  cali- 
bration points  for  each  instrument.  Since  the  position  indicator  is 
linear,  N  =  1  was  used  for  this  instrument.  The  thermistor  and  the 
conductivity  probe,  however,  required  N  =  2  and  N  =  4,  respectively 
(see  Figures  8  and  11).  The  polynomial  equations  were  then  used  within 
the  computer  program  to  convert  instrument  voltages  to  physical  values. 
In  the  conversion  of  conductivity  voltage  to  density,  corrections  were 
made  for  changes  in  temperature.  Figure  12  shows  the  effect  of  temper- 
ature on  the  conductivity  voltage  for  two  solutions  with  different 
salinities.  The  shape  of  the  hot  film  calibration  curves  (see  Fig- 
ure 13)  and  the  fact  that  below  a  critical  velocity  the  velocity  can 
no  longer  be  measured  accurately  make  it  desirable  to  fit  these  cali- 
bration points  with  a  sequence  of  straight  lines  rather  than  a  poly- 
nomial. If  the  hot  film  voltages  dropped  below  the  critical  values 
(points  1  in  Figure  13)  more  than  five  percent  of  the  time,  the  data 
were  discarded.  Variations  in  the  fluid  temperature  and  density  were 
accounted  for  when  the  velocity  was  computed. 

In  the  reduction  of  the  CST  data,  it  was  desirable  to  have  the 
data  output  at  nearly  uniform  intervals  of  about  0.5  mm  (vertical 
separation  distance).  To  accomplish  this,  m  values  of  sampled  data 
were  averaged  to  produce  one  output  point.  The  value  of  m  depended  on 
the  speed  of  the  traverse  at  the  time,  but  it  usually  fell  between 
10  and  30.  This  numerical  filter,  which  reduced  the  effects  of  high 
frequency  electrical  noise  (specifically  60-cycle),  was  applied  to  all 
three  channels  of  CST  data  (position,  temperature,  and  conductivity). 

Determination  of  the  Vaisala  frequency  at  a  point  requires  the 
density  gradient,  dp/dz,  at  the  point.  The  nature  of  differentiation 
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is  such  that  the  calculated  value  of  the  density  gradient  can  fluctuate 
greatly  due  to  small  errors  in  the  value  of  the  density  and/or  the 
position.  To  lessen  this  problem,  the  density  gradient  at  a  point  was 
calculated  by  first  fitting  a  least-squares  quadratic  to  five  points  in 
the  density  profile  (the  point  of  interest  and  two  points  on  .either 
side).  The  density  gradient  was  then  calculated  by  evaluating  the 
derivative  of  the  quadratic  at  the  point  of  interest.  The  thickness 
of  the  interface  6  was  determined  from. 


p  NMAX2 
where  NMAX  is  the  largest  value  of  the  Vaisala  frequency  taken  from 
the  core  of  the  interfacial  region. 

The  final  result  from  the  CST  data  was  a  tabulation  of  tempera- 
ture, specific  gravity,  and  Vaisala  frequency  versus  position.  The 
computer  program  then  plotted  these  profiles  using  the  Gould 
(series  5100)  plotter.  The  program  also  printed  out  the  position  of 
the  interface,  the  value  of  the  density  which  defines  the  interface, 
and  the  exact  time  that  the  sensors  crossed  the  interface. 

The  raw  data  from  the  ST  runs  consisted  of  four  or  five  time 
series:  thermistor  voltage,  conductivity  voltage,  and  the  output 
voltage  from  two  or  three  hot  film  anemometers  versus  time.  The  first 
step  in  reducing  these  data  was  to  convert  each  value  of  thermistor 
voltage  to  temperature  in  degrees  Fahrenheit.  The  temperature  at  each 
point  in  time  was  then  used  to  convert  the  conductivity  voltage  to  the 
fluid  density  (the  calibration  polynomials  discussed  above  were  used 
to  make  both  these  conversions).  Transformation  of  the  hot  film  volt- 
ages to  instantaneous  velocities  requires  knowledge  of  the  instantaneous 
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fluid  temperature  and  density  around  the  hot  film  sensors  (see 
Figure  13).  It  should  be  noted  that  the  thermistor,  the  conductivity 
probe,  and  the  hot  film  probes  were  separated  horizontally  by  a  small 
but  finite  distance  (  =  2  cm).  This  horizontal  separation  meant  that 
high  frequency  fluctuations  in  the  temperature  and  conductivity  at  the 
respective  probe  positions  were  not  related  in  general  to  the  high 
frequency  fluctuations  in  temperature  and  density  at  the  hot  film 
sensors.  If  these  fluctuations  are  not  filtered  out  they  may  induce 
fluctuations  in  the  velocity  measurements  which  do  not  actually  exist. 
Also,  since  the  temperature  and  density  corrections  were  applied  uni- 
formly to  both  channels  of  the  hot  film  data,  these  fluctuations  may 
have  resulted  in  unrealistic  values  of  the  velocity  correlations.  To 
reduce  the  effects  from  this  source  of  error,  the  temperature  and 
density  were  averaged  over  a  one  second  period  centered  about  the  time 
point  of  interest.  The  averaged  values  were  then  used  to  compute  the 
velocity  at  the  time  of  interest.  Mathematically,  this  process  is  as 
follows: 

1   i+n  -    1   i+n 

T.  =  i:    T(t=t.)  ,  p.  =  z    p(t=t.) 

2n+l  j=i-n     ^         2n+l  j=i-n     ^ 

V(t=t.)  =  V[HFV(t=t.),  T.,  p.]   . 

Where  HFV  is  the  hot  film  voltage,  V  the  velocity,  T  the  tempera- 
ture, and  p  the  fluid  density.  The  functional  relationship  between  V, 
HFV,  T,  and  p  is  shown  in  Figure  13  and  (2n+l)At  -   1  second. 

In  order  to  estimate  the  magnitude  of  this  potential  source  of 
error  a  simple  experiment  was  conducted.  Identical  time  series  records, 
for  a  point  in  the  interfacial  region,  were  reduced  by  two  different 
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methods.  With  the  first  method  the  instantaneous  velocity  and  all  the 
correlations  were  computed  using  a  constant  value  of  temperature  and 
density  (the  average  values  of  temperature  and  density  for  the 
540  second  time  period).  The  second  method  used  the  instantaneous 
values  of  temperature  and  density  to  compute  the  instantaneous 
velocities  and  the  velocity  correlations.  A  comparison  of  the  results 
showed  that  the  relative  difference  between  the  values  was  always  less 
than  0.1  percent.  It  was  concluded,  therefore,  that  fluctuations  in 
the  temperature  and  density  do  not  result  in  a  significant  source  of 
error. 
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CHAPTER  IV 
RESULTS 

Introduction 

The  data  presented  in  this  chapter  are  the  results  of  those 
experiments  considered  accurate,  consistent,  and  reliable.  The  pri- 
mary factor  in  making  this  determination  was  the  mean  velocity.  If 
the  mean  velocity  at  a  point  varied  significantly  over  a  long  period 
of  time  or  if  the  velocity  profile  away  from  the  interface  varied 
radically  with  position,  the  data  were  discarded  as  invalid. 

Figure  14  gives  a  chronological  listing  of  the  experiments  and 
the  corresponding  global  parameters  which  define  the  state  of  flow. 
Figure  15  lists  some  of  the  experimentally  derived  quantities  for  each 
test. 

Five  experiments  were  conducted  at  a  location  800  cm  downstream 
from  the  entrance  while  four  others  were  conducted  at  1600  cm.  The 
overall  density  difference  was  maintained  between  2.2  and  2,7  percent 
while  the  mean  velocity  varied  between  2.7  cm/sec  and  7.8  cm/sec.  The 

Reynolds  number  for  the  upper  layer,  based  on  the  equivalent  pipe 

4 
diameter,  was  always  greater  than  2.0  x  10  .  Thus,  the  upper  layer 

was  always  above  the  critical  Reynolds  number  for  turbulent  flow, 

while  the  flow  deep  within  the  bottom  layer  was  always  laminar. 
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Density  Structure 

The  physical  appearance  of  the  density  interface  is  similar  to 
the  surface  of  water  at  the  point  when  wind  first  starts  to  generate 
waves,  except  that  the  time  scales  are  longer.  Small  amplitude  wave- 
lets having  large  wave  numbers  and  short  crests  move  along  the  inter- 
face with  an  apparently  random  distribution  in  both  space  and  time. 
The  wave  lengths,  crest  lengths  and  amplitudes  of  these  waves  increase 
very  slowly  with  distance  downstream.  Many  of  the  wavelets  form  sharp 
crests,  break  and  eject  salt  water  into  the  upper  layer.  At  times, 
the  ejected  salt  water  falls  back  to  the  interface  before  diffusing 
completely.  The  structure  of  the  density  interface  and  the  dynamics 
of  the  mixing  process  clearly  indicate  that  this  flow  is  always  within 
regime  two  (see  Chapter  I).  Although  the  interface  is  continuously 
moving,  it  always  remains  distinct  and  easily  identifiable  with  no 
large  coherent  overturning  structures  present. 

Appendix  B  contains  the  quasi-instantaneous  profiles  of  density, 
temperature  and  Vaisala  frequency  obtained  by  the  CST  technique. 
Because  the  density  interface  is  in  a  state  of  agitation  there  is  con- 
siderable variance  between  density  profiles,  even  between  consecutive 
profiles  taken  under  the  same  conditions.  Some  of  the  profiles 
exhibit  what  appear  to  be  double  interfaces  or  unstable  regions. 
These  structures  may  be  aberrations  resulting  from  the  passage  of  a 
wave  during  the  10  to  30  seconds  required  for  the  probes  to  traverse 
the  interfacial  region  or  they  may  actually  exist.  When  the  structures 
exist,  they  appear  in  all  profiles  taken  during  the  particular  experi- 
ment. When  the  structures  result  from  the  passage  of  a  wave,  they 
appear  and  disappear  at  random.  The  density  plots  from  October  26,  1978, 
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and  May  3,  1977,  illustrate  the  first  case  while  those  from 
October  13,  1978,  illustrate  the  latter.  The  small  perturbations  in 
the  density  profiles  found  deep  within  each  layer  are  most  likely  the 
result  of  noise  on  the  conductivity  signal.  These  fluctuations,  most 
pronounced  in  the  salt  layer  where  the  conductivity  probe  is  least 
sensitive,  translate  into  large  fluctuations  in  the  Vaisala  frequency 
profiles. 

Due  to  the  large  variance  between  individual  profiles,  a  more 
meaningful  measure  of  the  density  structure  is  obtained  by  first  non- 
dimensionalizing  the  density  and  Vaisala  frequency  profiles,  and  then 
forming  an  ensemble  average  by  overlaying  several  of  these  plots. 
Figures  16,  18,  20,  22  and  24  show  the  overlayed  density  profiles 
while  Figures  17,  19,  21,  23  and  25  show  the  corresponding  Vaisala 
frequency  plots.  The  nondimensional  density,  p*,  is  defined  by 

p(s/6)  -  p„ 


P*U/&) 


^^i  -  %^ 


where  p  is  the  density  of  the  upper  layer  and  p  the  density  of  the 
lower  layer.  The  vertical  coordinate  5  is  nondimensionalized  by  the 
interfacial  thickness  parameter,  6  defined  by 

2g  Ap 


=  (p^  +  pj   NMAX^ 

where  g  is  the  acceleration  of  gravity  and  NMAX  is  the  largest  value 
of  Vaisala  frequency  from  the  core  of  the  interfacial  region.  The 
dimensionless  Vaisala  frequency,  N*,  is  obtained  by  dividing  the 
Vaisala  frequency  by  NMAX.  The  solid  line  on  each  of  the  nondimen- 
sionalized density  plots  is  given  by 
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9*U/&)   =  -^  M  -  tanh  (^n 


while  the  solid  line  on  the  nondimensionalized  Vaisala  frequency  plots 
is  given  by 

N*(5/<S)  =  sech  (2^/6) 
These  functional  representations  of  the  density  and  Vaisala  frequency 
are  generally  a  good  fit  to  the  experimental  data.  The  greatest 
deviation  usually  occurs  in  the  region  above  the  interface  where  the 
curvature  of  the  density  profile  is  maximum. 

Figure  26  shows  a  plot  of  the  interfacial  thickness  parameter,  6, 
versus  the  Keulegan  number,  K.  The  solid  dots  represent  the  data 
taken  at  Station  800  while  the  open  circles  are  data  taken  at  Sta- 
tion 1600.  The  solid  line  is  the  line  of  least  square  error  through 
the  solid  dots.  The  corresponding  correlation  coefficient  is  0.43, 
which  for  25  degrees  of  freedom,  is  significant  at  the  0.05  level. 
The  slope  of  this  line  gives 

^  =  0.0126  cm. 

The  data  sets  defined  by  0£K<4  indicate  a  positive  correlation  also 
exists  between  &   and  x.  The  average  interfacial  thickness  for  the 
data  in  this  range  was  0.44  at  Station  800  while  at  Station  1600  it 
was  0.64,  thus 

^=  2.5  X  10""^ 
Ax 

The  correlation  coefficient  between  6  and  x  for  these  data  sets  was 
0.59,  which  given  the  28  degrees  of  freedom  is  again  significant  at 
0.001  level.  Due  to  the  limited  number  of  stations  and  experiments 
it  was  not  possible  to  determine  the  functional  relationship  between 
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6  and  X  or  between  6  and  K,  but  there  is  a  positive  correlation  in 
both  cases. 

Mean  Velocity  *'.. 

The  mean  velocity  profiles  fall  into  two  distinctly  different 
groups.  The  primary  feature  differentiating  the  groups  in  the  bound- 
ary layer  thickness,  6  ,  defined  as  the  point  where  the  mean  velocity 
reaches  95  percent  of  the  freestream  velocity.  The  freestream  velocity, 
U,  is,  in  turn,  defined  as  the  average  velocity  15  cm  or  more  above  the 
density  interface.  The  velocity  in  this  region  is  \jery   nearly  uniform. 
Figure  27  shows  the  nondimensionalized  velocity  profiles  for  the  first 
group;  the  boundary  layer  thickness  for  these  profiles  is  between 
10.5  cm  and  11.5  cm.  In  the  second  group.  Figure  28,  the  boundary  layer 
thickness  is  only  between  1.15  cm  and  3.2  cm.  The  two  nondimensional- 
ized profiles  also  have  different  shapes.  For  comparison,  the  two 
velocity  profiles  shown  in  Figures  27  and  28  are  drawn  together  in 
Figure  29.  In  the  region  (5-n)/6  >  0.6  the  two  profiles  are  identical, 
within  experimental  error.  However,  below  this  region  the  profiles  are 
distinctly  different. 

No  statistically  significant  distinction  can  be  made  between  pro- 
file groups  on  the  basis  of  any  of  the  parameters  listed  in  Figure  14, 
However,  there  is  a  slight  tendency  for  the  group  two  profiles  to  occur 
at  lower  velocities  and  at  locations  closer  to  the  entrance.  The  small 
variations  in  the  boundary  layer  thickness  for  the  profiles  in  Figure  27 
have  no  significant  correlation  with  mean  velocity,  distance  downstream 
or  the  Reynolds  number  based  on  these  two  quantities.  Therefore,  the 
flow  is  believed  to  be  fully  developed  and  in  equilibrium.  This, 
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however,  is  not  the  case  for  the  group  two  profiles  (Figure  28).  The 
boundary  layer  thickness  for  this  set  correlates  well  with  the  inverse 
of  the  Reynolds  number  based  on  the  downstream  position.  The  correla- 
tion coefficient  is  0.969  which,  given  two  degrees  of  freedom,  is  still 
significant  at  the  0.03  level. 

The  origin  of  the  vertical  coordinate  in  both  Figures  27  and  28  is 
defined  as  the  point  where  the  velocity  equals  50  percent  of  the  free- 
stream  velocity.  This  origin  is  shifted  a  small  amount,  n,  from  the 
origin  defined  by  the  density  interface  (i.e.  5=0),  This  small  shift 
between  the  density  and  velocity  profiles  appears  reasonable  upon  exam- 
ination of  the  individual  density  profiles  (Appendix  B).  A  small  dis- 
crepancy of  less  than  0.2  cm  could  result  from  errors  in  positioning 
the  conductivity  and  hot  film  probes,  while  larger  shifts  are  usually 
associated  with  the  existence  of  a  double  interfacial  structure.  It 
appears  that  the  velocity  profile,  in  some  way,  "locks  in"  on  the  upper 
most  density  interface. 

The  experiments  which  exhibit  the  group  two  velocity  profiles  also 
exhibit  relatively  smooth  density  profiles  having  a  single  density 
interface.  On  the  other  hand,  the  experiments  which  fall  into  group  one 
have  density  profiles  which  exhibit  either  a  double  interface  or  a 
highly  agitated  structure  (see  Appendix  B).  This  correlation  between 
velocity  profiles  and  density  structure  extends  to  the  vertical  shift 
parameter,  n-  As  shown  in  Figures  27  and  28  the  group  one  velocity  pro- 
files have  values  of  n  ranging  from  0.38  to  0.80  cm  while  for  the  group 
two  profiles,  n  falls  between  0.05  cm  and  0.25  cm. 

The  mean  velocity  profile  for  group  one  is  plotted  on  log-normal 
paper  in  Figure  30.  This  illustration  points  out  the  existence  of  a 
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region  above  the  interface,  (?  -  n)/&^     >   0.15,  where  the  velocity  pro- 
file has  a  logarithmic  shape.  The  group  two  profiles  have  no  such 
region. 


The  Richardson  Number 

The  gradient  Richardson  number  at  any  point  is  defined  by: 


Ri 


I    3U  x^ 


N2 


■(1E7  ' 


where  N  is  the  Brunt-Vaisala  frequency.  Given  the  Vaisala  frequency 
profile. 


where  by  definition. 


N(s)  =  NMAX  sech(|^)  , 


NMAX2  =  ^SlAp. 


^  a    u' 


u  =  u  f(  p)   , 

u 


and  the  velocity  profile, 

Uf(p 

U 

The  Richardson  number  profile  can  be  expressed  as: 

2g  Ap  62 


Ri"  (0 


(P/P,)U26 


sech  (|^) 


u 


where  f'(-^^)  is  the  derivative  of  f  with  respect  to  ^^  . 
u  u 

The  gradient  Richardson  number  at  the  inflection  point  in  the  velocity 

profile  (?=n)  then  becomes: 
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Rig(5=r,) 


2g  Ap  62 
U2(P,+P,)6 


sech(|l) 
f  (0) 


The  gradient  Richardson  number  at  5=n  was  calculated  for  each  experi- 
ment using  the  above  equation  and  the  appropriate  values  listed  in 
Figures  14  and  15  (the  average  value  of  the  6's  listed  in  Figure  15 
was  used  for  6).  The  values  of  f'(0)  were  obtained  for  the  two 
velocity  groups  from  Figures  27  and  28.  The  results  are  shown  below: 


Exp.  No. 


Date 


Group 


Ri'  (C=n) 


f  (0) 


5/3/77 

M 

2 

0.08 

8.0 

3 

10/19/77 

0.13 

8.0 

6 

10/26/78 

0.52 

8.0 

7 

10/26/78 

1.06 

8.0 

8 

11/8/78 

0.07 

8.0 

1 

7/6/76 

2 

88.60 

0.8 

4 

7/24/78 

2 

1.23 

0.8 

5 

10/13/78 

2 

7.18 

0.8 

9 

12/6/78 

2 

7.06 

0.8 

Thus 


and 


Ri  (C=n)  1  1.06  for  group  one  profiles 


Ri  (C=n)  1  1.23  for  group  two  profiles. 


Shear  Stress  Profiles 

The  total  shear  stress  in  a  turbulent  flow  is  the  sum  of  the 
viscous  and  the  Reynolds  stresses  defined  as 

8U 


■'T  =  \  ^  H  "  ^  W  ''"'^'   ' 


where  Tx  is  the  total  stress,  x  is  the  viscous  stress,  and  x. 
T  V  t 

is  the  Reynolds  stress. 


56 


Figures  31  and  32  show  the  total  shear  stress  profile  correspond- 
ing to  the  group  one  and  group  two  velocity  profiles  respectively.  The 
stress  has  been  normalized  by  x-,  the  total  stress  at  the  velocity 
interface,  5=n.  Both  these  profiles  exhibit  a  marked  increase  in  the 
neighborhood  of  the  interface,  dropping  off  rapidly  on  either  side. 
The  fact  that  the  total  shear  stress  beomes  very  small  just  below  the 
density  interface  is  in  agreement  with  the  observation  of  little  or  no 
motion  in  the  bottom  layer  even  after  the  upper  layer  has  been  moving 
for  three  to  six  hours.  The  rapid  decrease  in  the  stress  profile  above 
the  interface,  however,  is  unexpected  since  equations  2.2  and  2.4  pre- 
dict a  linear  stress  profile  between  the  freesurface  and  the  density 
interface.  In  an  attempt  to  explain  this  unexpected  behavior,  the  theo- 
retical shear  stress  profile  was  calculated  from  equations  2.1  through 
2.3.  These  equations  account  for  the  additional  effects  of  side  wall 
friction,  mean  advective  accelerations  and  mean  density  gradients.  The 
latter  two  result  from  boundary  layer  growth  and  turbulent  diffusion  at 
the  interface,  set-up  of  the  density  interface  and  set-down  of  the 
freesurface.  Details  of  the  solution  techniques  are  presented  in 
Appendix  E.  Figure  33  shows  the  computed  shear  stress  profiles  under 
group  one  conditions  and  the  magnitude  of  the  five  parameters  which 
determine  the  shape  of  the  stress  profiles.  These  profiles  more  closely 
represent  the  data  (see  Figure  31)  than  does  the  simple  linear  model. 
However,  this  solution  is  still  not  completely  satisfactory.  First, 
the  calculated  stress  is  only  40  to  50  percent  of  the  measured  stress 
at  the  interface.  In  addition,  except  for  profiles  one  and  two, 
the  magnitude  of  the  parameters  are  unreasonably  large,  A  possible 
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explanation  for  why  this  analysis  is  unable  to  predict  the  shear  stress 
profile  is  given  in  Chapter  V. 

The  interfacial  friction  factor,  defined  by 

8  T, 

f. 


^   pU2 

is  plotted  versus  Reynolds  number  in  Figure  34  along  with  the  line 
corresponding  to  turbulent  flow  over  a  hydraulically  smooth  rigid 
surface.  This  figure  shows  that,  at  the  same  Reynolds  number,  the  shear 
stress  at  the  interface  is  generally  less  than  or  equal  to  the  stress 
at  a  smooth  rigid  surface.  The  data  points,  except  for  one,  cluster 
around  a  line  with  a  -2  slope  which  implies  the  interfacial  shear  stress 
is  independent  of  the  mean  flow  velocity.  No  significant  correlation  is 
found  between  the  interfacial  friction  factor  and  the  data  groups. 

Figures  35  and  36  show  the  ratio  of  viscous  stress  to  total  stress 
for  the  group  one  and  two  velocity  profiles.  There  is  again  an  obvious 
difference  between  groups.  For  group  one,  the  viscous  stress  drops  off 
in  an  exponential  manner  above  the  interface;  at  the  interface,  the 
viscous  stress  is  only  60  to  70  percent  of  the  total  stress.  For  group 
two,  the  viscous  stress  remains  large  above  the  interface  until 
(5-n/6  )  >  1.0  and  then  drops  off  rapidly.  For  (c-n/^  )  <  1.0  the 
viscous  stress  accounts  for  nearly  the  entire  stress  and  in  some  cases 
balances  negative  Reynolds  stresses.  Again  for  the  purpose  of  compari- 
son, the  curves  from  Figures  35  and  36  are  shown  together  in  Figure  37. 

Energy  Spectra 

Energy  spectra  of  the  horizontal  and  vertical  velocity  fluctua- 
tions were  measured  during  four  of  the  nine  experiments  at  several 
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positions  above  the  interface.  These  plots,  along  with  pertinent  facts 
on  the  measurement  technique,  are  given  in  Appendix  D.  The  spectra  mea- 
sured well  above  the  interface  generally  have  a  smoothly  varying  shape 
indicating  broad  band  phenomena  denoting  a  turbulent  flow.  The  energy 
content  of  the  high  frequency  fluctuations  appears  to  be  about  the 
same  for  the  horizontal  and  vertical  components  indicating  local 
isotropy.  However,  at  lower  frequencies  the  energy  in  the  vertical 
fluctuations  drops  below  that  for  the  horizontal  fluctuations  due  to 
the  presence  of  the  stable  density  interface. 

The  spectra  measured  close  to  or  in  the  density  interface  exhibit 
a  more  erratic  behavior  and  a  rather  sudden  change  in  the  slope  of  the 
energy  curve  between  one  and  two  Hertz;  this  corresponds  roughly  to 
the  maximum  Vaisala  frequency  of  the  interface.  There  is  often  a  marked 
rise  in  the  energy  level  in  this  frequency  band,  especially  in  the  ver- 
tical velocity  component  and  a  rapid  decrease  in  the  energy  at  higher 
frequencies. 

This  characteristic  of  the  energy  spectra  indicates  the  presence 
of  internal  wave  energy.  The  fact  that  the  spectra  drops  off  rapidly 
above  the  Vaisala  frequency  is  indicative  of  the  fact  that  disturbances 
with  frequencies  above  the  Vaisala  frequency  are  damped  rapidly.  The 
density  interface  prevents  the  transfer  of  high  frequency  energy  from 
one  side  of  the  interface  to  the  other. 

Finally,  the  spectra  measured  on  December  6,  1978,  at  c  =  0.37  cm 
and  on  July  6,  1976,  at  5  =  0.46  cm  correspond  to  group  two  velocity 
profiles.  These  spectra  indicate  the  presence  of  internal  waves. 


59 


Ln 

U3 

■=:*• 

CM 

O 

00 

>* 

o 

«=f 

? 

CO  en 

■=J- 

<- 

CO 

r^ 

r^ 

00 

UD 

Ln 

CO 

Q. 

\=>  o 

• 

• 

■ 

• 

• 

• 

• 

• 

> 

< 

Q.I — 

CVJ 

lO 

CM 

o 

CM 

o 

CM 

o 

CO 

OT 

X 

f— 

QJ 

CM 

O 

<JD 

r^ 

1^ 

CO 

cn 

Ln 

Q 

^ 

• 

• 

• 

• 

■ 

• 

• 

Q. 

1=0 

CO 

<y\ 

UD 

r— 

00 

LO 

r-- 

^— 

^— 

< 

Q. 

r^ 

<JD 

O 

•vl- 

en 

«=!- 

00 

CO 

CO 

cn 

CSJ 

r— 

CO 

O 

^o 

O 

o 

o 

CM 

CO 

1 — 

o 

1 

a 

• 

■ 

• 

• 

• 

• 

• 

• 

X 

o 

<X> 

r~ 

"=1- 

CO 

o 

CO 

"^ 

lO 

V£l 

p  1— 

Lf) 

lO 

<X) 

U3 

«=1- 

CM 

^ 

Ln 

1=3 

CM 

CM 

ro 

LO 

CO 

o 

r~~ 

CM 

LO 

X 

«^ 

f> 

"^""^ 

' 

' 

\S^ 

CM 

"^ 

CO 

LO 

Ln 

en 

r^ 

Ln 

<U 

<d- 

r^ 

CM 

o 

t~^ 

en 

CO 

r^ 

CvJ 

Q 

■> 

cy> 

«s 

ro 

Ln 

en 

CO 

CO 

1=3 

^— 

U3 

CO 

r~ 

<JD 

CM 

00 

CO 

CO 

CM 

CM 

ro 

Ln 

CO 

Ln 

CO 

<X) 

CO 

5^ 

o 

O 

Ln 

^ 

o 

o 

o 

o 

LD 

CM 

CM 

•JD 

Ln 

Ln 

Ln 

^ 

r~~ 

Q.       3 

• 

• 

• 

• 

• 

• 

• 

• 

• 

O     Q. 

OJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o 
ai 
on 

CM 

(_) 

o 

CM 

O 

o 

O 

CM 

<^ 

en 

Ln 

E 

1=)   B 

a; 

r-. 

CM 

1— 

Ln 

Ln 

r^ 

1^ 

o 

u 

u 

(/) 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CM 

CO 

"^ 

1 — 

•* 

U3 

"=1- 

r^ 

'^f 

o 
o 
II 

X  e 

O 

o 

o 

CD 

o 

O 

o 

o 

o 

■> 

o 

O 

o 

o 

O 

o 

O 

o 

o 

o 

00 

00 

00 

00 

00 

VD 

<X3 

LD 

U3 

o 

ro 

>* 

,-1 

VD 

o^ 

CM 

CO 

UD 

yD 

CO 

LU 

V£) 

CO  I~-. 

r—  r^ 

00 

n-  CO 

CM  CO 

CM  00 

00  CO 

<J3  CO 

1 — 

>,r^ 

r-. 

r^ 

>,r-. 

r^ 

r- 

r^ 

r~ 

r^ 

II 

cC 

1—  CTt 

>,CT. 

4->  a> 

r-^O^ 

+->  en 

+->  (Tl 

4->  en 

>  en 

O  C7> 

Q 

3  1 — 

n3  1 — 

O  1 — 

3   r — 

U  1— 

o  <— 

O  r— 

O  1 — 

O)  1 — 

X 

■-3 

s: 

o 

'n 

o 

o 

O 

■zi 

Q 

en 
II 

. 

Q.      • 

X    O 

, — 

CM 

ro 

^ 

Ln 

LD 

1 — 

CO 

CTi 

a 

UJ  ^ 

Q 

00 

en 


a. 


00 

o 
—1 

CD 
Q 


CC 


Q. 
X 


CD 


60 


o 

en 

^ 

LO 

CO 

O 

d 

o 

^-— «K 

(J 

CO 

CO 

o 

o 

CO 

>* 

o 

CO 

CO 

LO 

Lf) 

LO 

LO 

«^ 

o 

«^ 

1— 

o 

• 

o 

• 

o 

d 

<~> 

d 

• 

• 

o 

<_> 

-^ 

n 

t\ 

M 

«\ 

9\ 

»> 

r* 

w 

OJ   B 

OJ  f— 

r^ 

o 

LO  CM 

■v^  CO 

en 

00 

00 

^J- 

S-    o 

^  LO 

«d- 

«D 

CO  >* 

CO  00 

r^ 

CO 

vo 

^ 

3  s_ 

,       . 

• 

• 

•              • 

•         • 

• 

• 

• 

• 

(/)  »4- 

o  o 

o 

O 

o  o 

o  o 

o 

o 

o 

o 

ro 

«t 

«^ 

r—   CO 

CO 

, 

CO 

CO  «d- 

o 

vo 

r^ 

vo  r^ 

<£i 

^ 

CO 

LO  «^ 

r-~ 

r» 

Lf) 

"^ 

a            • 

, 

a 

a 

•         • 

• 

• 

• 

• 

o  o 

o 

o 

o 

o  o 

o 

o 

o 

o 

•* 

«( 

•» 

#> 

« 

#t 

ro 

Od 

Cvl 

00 

r^ 

vo 

VO 

CO 

LO 

vo 

vo 

"=3- 

O 

o 

d 

d 

o 

o 

■> 

ro* 

o 

\n 

o 

U3 

^1- 

o 

CM 

^ 

CO 

Q. 

3 

ir> 

<ys 

r^ 

C\J 

^— 

CO 

r^ 

00 

CO 

< 

• 

• 

• 

• 

• 

• 

• 

• 

• 

ct[ 

Q. 

,^ 

U3 

CTl 

00 

r^ 

f— 

en 

<n 

1^ 

1 — 

OJ 

i — 

1 — 

1 — 

CM 

CM 

' — 

CO 

'~ 

o 

O 

CVJ 

in 

in 

00 

CO 

LO 

LO 

cr  E 

CO 

vo 

C\J 

n— 

CO 

CO 

VO 

o 

o 

• 

a 

a 

a 

• 

• 

• 

• 

• 

o 

O 

O 

CD 

o 

o 

o 

o 

o 

o 

LO 

o 

in 

CVJ 

o 

00 

«* 

vo 

CO 

LO 

lO 

r— 

3   E 

• 

• 

• 

• 

• 

• 

• 

• 

• 

<<3         (J 

OO 

^ 

o 

o 

o 

^ 

en 

>=r 

O^ 

CO 

o-i 

<:*• 

r~» 

00 

CO 

CVJ 

o 

CM 

o 

r— 

o 

r^ 

■  r- 

o 

o 

o 

o 

o 

o 

o 

o 

o 

M- 

, 

a 

a 

a 

a 

a 

• 

• 

• 

o 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

UD 

CM 

CO 

«* 

CM 

CO 

CM 

o 

CO 

o 

CJ 

^3- 

^ 

CM 

o 

CM 

CT> 

*E    O) 

OsJ 

CvJ 

CM 

CM 

CM 

CM 

CM 

f— 

3    O    to 

• 

• 

« 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

O 

00 

LO 

Q.             (N 

r^ 

o 

f_ 

00 

en 

o 

r— 

CM 

r^ 

--"^ 

"--^     CM 

o 

<d- 

IT) 

ir> 

LO 

LO 

^ 

LO 

CO 

G 

•r-     E 

<u 

o 

o 

o 

o 

o 

o 

o 

o 

o 

'^^'N 

h-         O 

CO 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•i 

o 

o 

o 

o 

o 

o 

o 

C3 

o 

tl 

, 

Q.     • 

X    O 

^— 

CvJ 

CO 

«* 

LO 

vo 

t^ 

OO 

en 

• 

UJ  Z 

M- 

CO 


ID 

cy 

Q 
LU 

OL 
ID 
OO 

<: 


Q 

OO 


Q. 
X 


a: 


61 


DATE:     IP/fi    /7R 
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FIGURE  16:  NONDIMENSIONAL  DENSITY  PROFILE  12/6/78 
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CHAPTER  V 
CONCLUSIONS 

This  study  has  documented  the  structure  of  the  velocity  and 
density  transition  region  in  a  two-layer  flow  by  experimentally  mea- 
suring the  time  mean  velocity  and  density  profiles,  the  turbulent  and 
viscous  shear  stress  and  the  turbulent  energy  spectra.  In  addition, 
the  measured  data  were  used  to  calculate  the  interfacial  friction  factor 
and  the  gradient  Richardson  number  at  the  inflection  point  of  the  veloc- 
ity profile.  From  these  measurements,  the  following  conclusions  are 
reached. 

1)  In  a  weakly  turbulent  stratified  shear  flow,  the  interfacial 
region  exists  in  one  of  two  measurably  different  states.  This  bimodal 
character  is  not  manifest  in  the  magnitude  of  the  flow  properties,  but 
in  the  distribution  of  these  properties  through  the  interfacial  region, 

2)  No  parameter  or  dimensionless  number  derived  from  global  or 
overall  flow  properties  accurately  predicts  the  existence  of  the 
two  modes  or  groups.  However,  the  gradient  Richardson  number  at  the 


The  term  "global"  in  this  context  implies  a  integral  scale  property, 
for  example  the  mean  freestream  velocity,  the  depth  or  width  of  the 
channel,  the  total  density  jump  across  the  density  interface,  or  the 
Richardson  number  based  on  the  above.  This  is  in  contrast  to  the  term 
"local"  which  refers  to  point  properties  such  as  the  density  or  the 
gradient  Richardson  number  based  on  the  velocity  and  density  gradient 
at  a  point  within  the  flow. 
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inflection  point  of  the  velocity  profile  is  distinctly  different  for 
each  group.  In  group  one  the  Richardson  number  is  less  than  or  equal  to 
one,  while  in  group  two  the  converse  is  true. 

3)  The  turbulent  energy  spectra  measured  within  the  boundary  layer 
above  the  density  interface  indicate  the  presence  of  internal  wave 
motions  under  group  two  conditions.  Under  group  one  conditions,  no  in- 
ternal wave  subrange  is  indicated  by  the  energy  spectra. 

4)  The  total  shear  stress  in  the  x-y  plane  is  small  everywhere  ex- 
cept in  a  narrow  region  around  the  density  interface  where  the  stress 
reaches  a  maximum  (Figures  31  and  32). 

5)  Over  the  range  of  conditions  investigated,  the  interfacial  fric- 
tion factor  is  generally  less  than  or  equal  to  the  friction  factor  for 

a  hydraulically  smooth  surface  at  the  same  Reynolds  number.  Furthermore, 
the  friction  factor  is  proportional  to  Re   (see  Figure  34).  Therefore, 
the  interfacial  shear  stress  is  independent  of  the  freestream  velocity 
when  the  Reynolds  number  is  between  30,000  and  65,000. 

The  bimodal  character  of  the  flow  near  the  interface  is  indicated 
by  the  mean  velocity  and  shear  stress  profiles,  the  profiles  of  viscous/ 
total  shear  stress,  and  the  gradient  Richardson  number  at  the  inflection 
point  of  the  velocity  profile.  For  each  of  the  above,  the  experimental 
data  divide  into  one  of  two  distinct  groups. 

With  respect  to  the  mean  velocity,  there  are  two  distinguishing 
characteristics.  First,  the  group  one  profiles  have  a  boundary  layer 
thickness  on  the  order  of  10  cm  while  group  two  profiles  have  a  thick- 
ness of  only  one  to  three  centimeters.  Secondly,  when  nondimensional- 
ized,  the  velocity  profiles  for  the  two  groups  have  significantly 
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different  shapes.  The  group  one  profiles  have  a  shape  characteristic 
of  a  turbulent  boundary  layer  over  a  rigid  wall,  possessing  a  log-linear 
distribution  over  a  limited  region.  The  group  two  profiles  do  not  ex- 
hibit this  characteristic  shape.  The  nondimensionalized  shear  stress 
profiles  also  bifurcate,  each  group  having  a  distinctly  different  shape. 

The  profiles  of  viscous/total  shear  stress  point  again  to  the 
bimodal  character.  For  group  one,  the  stress  is  35  to  40  percent  turbu- 
lent at  the  interface  and  the  viscous  component  drops  rapidly  to  less 
than  10  percent  before  reaching  the  outer  edge  of  the  boundary  layer 
(see  Figure  37).  The  group  two  data,  however,  show  the  shear  stress  to 
be  about  90  percent  viscous  throughout  the  boundary  layer  falling  off 
only  outside  this  region.  In  addition,  the  data  for  group  one  have 
little  variance  while  group  two  data  are  widely  scattered. 

The  density  profiles  also  show  a  character  that  correlates  with  the 
velocity  groups.  The  group  one  experiments  almost  invariably  exhibit  an 
agitated  density  interface,  many  times  possessing  two  or  three  distinct 
interfaces  while  the  group  two  experiments  usually  have  a  smoothly  chang- 
ing density  profile. 

The  fact  that  no  global  parameter  or  dimensionless  number  correlates 
significantly  with  the  two  groups  indicates  that  the  existence  of  the 
two  dynamic  states  is  not  determined  by  the  global  properties  of  the 
flow.  Even  the  interfacial  friction  factor,  which  is  a  measure  of  the 
maximum  shear  stress  at  the  interface,  shows  no  significant  correlation 
with  the  groups.  In  contrast,  the  gradient  Richardson  number  at  the 
inflection  point  of  the  velocity  profile  is  found  to  have  a  strong 
correlation  with  the  two  groups.  This  verifies  the  fact  that  the  mean 
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velocity  and  density  distributions  through  the  transition  region  are 
significantly  different  from  one  group  to  the  other.  Since  the  gradient 
Richardson  number  is  calculated  from  the  velocity  and  density  distribu- 
tion, it  can  not  be  considered  as  the  factor  that  dictates  which  dynamic 
state  will  exist,  but  must  be  viewed  only  as  an  indicator  of  which  state 
exists  at  the  time.  The  significance  of  the  gradient  Richardson  number 
at  the  inflection  point  of  the  velocity  profile  is  brought  into  focus 
when  combined  with  two  other  observations.  First,  the  velocity  inflec- 
tion point  is  always  at  or  above  the  density  inflection  point.  Second, 
the  density  transition  region  is  always  thinner  than  the  velocity  tran- 
sition region  (at  least  for  this  flow  regime).  These  two  factors  com- 
bined imply  that  the  gradient  Richardson  number  everywhere  above  the 
velocity  inflection  point  is  less  than  or  equal  to  the  Richardson 
number  at  the  velocity  inflection  point.  Therefore,  the  gradient 
Richardson  number  is  less  than  or  equal  to  one  everywhere  above  the  in- 
flection point  for  group  one.  For  group  two,  this  is  not  the  case.  To 
maintain  turbulence,  the  gradient  Richardson  number  must  be  less  than 
or  equal  to  one  (see  Turner,  1973). 

The  author  suggests  the  above  facts  point  to  the  following  conclu- 
sions. First,  the  two  groups  represent,  in  one  case,  a  well  developed 
turbulent  boundary  layer  (group  one)  and,  in  the  second,  a  gravitation- 
ally  stabilized  boundary  layer  (group  two).  This  conclusion  is  sup- 
ported by  the  fact  that  the  turbulent  energy  spectra,  measured  within 
the  boundary  layer  above  the  density  interface,  indicate  the  existence 
of  an  internal  wave  subrange  for  group  two  experiments.  The  energy 
spectra  measured  in  the  same  region  (although  not  at  exactly  the  same 
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point)  show  no  such  subrange  for  group  one  experiments.  Second,  since 
no  correlation  is  found  between  the  dynamic  modes  (groups)  and  any  of 
the  global  flow  parameters,  the  author  believes  the  data  groupings 
must  be  related  to  the  initial  conditions  at  the  start  of  each  experi- 
ment. Although  there  was  no  conscious  or  observable  difference  in  the 
initiation  of  each  experiment,  subtle  differences  may  have  been  present. 

The  shear  stress  profiles  for  both  groups  (Figures  31  and  32)  have 
very   unexpected  shapes  which  are  difficult  to  explain.  The  stress  is 
small  everywhere  except  in  a  narrow  region  around  the  interface,  where 
the  magnitude  increases  rapidly.   One  to  three  centimeters  below  the 
interface  the  shear  stress  goes  to  zero.  This  is  in  agreement  with  the 
observation  of  little  or  no  motion  below  the  interface,  even  after  the 
flow  in  the  upper  layer  has  been  in  motion  for  three  to  six  hours. 

An  attempt  was  made  to  solve  equations  2.1  through  2.3  in  order  to 
determine  the  origin  of  the  unexpected  stress  profile.  The  equations 
were  solved  (see  Appendix  E)  under  the  assumption  that  the  velocity  and 
density  profiles  were  each  self-similar  and  that  each  required  only  one 
length  scale  for  nondimensionalization.  This  approach  failed  to  repro- 
duce accurately  the  measured  shear  stress  profiles  (compare  Figures  31 
and  33).  One  can  only  assume  that  significant  terms  are  omitted 
from  equations  2.1  through  2.3  (terms  which  were  assumed  small,  see 
Appendix  E)  or  that  the  self-similarity  assumption  is  incorrect. 
The  author  believes  the  latter  is  more  likely,  since  nondimensionalizing 
the  velocity  profile  by  one  length  scale  implies  the  boundary  layer 
above  the  interface  grows  at  the  same  rate  as  the  boundary  layer  below 
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the  interface.  This  may  not  be  the  case,  because  the  flow  above  the 
interface  appears  turbulent  while  the  flow  below  the  interface  is 
laminar. 
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APPENDIX  A 
PROFILE  SHAPING  ELEMENT 

The  profile  shaping  element  (PSE)  was  designed  to  produce,  at  the 
entrance  to  the  test  section,  a  velocity  profile  which  closely  approxi- 
mates fully  developed  turbulent  flow  in  a  rectangular  channel,  thereby 
minimizing  longitudinal  variations  within  the  test  section.  The  fully 
developed  profile  is  closely  represented  by: 

Where    U^^^  =  U(y=0,  z=D), 

w  is  the  half-width  of  the  tank,  D  the  depth  of  the  channel,  y  the 
lateral  distance  from  the  center  line  of  the  tank,  and  z  the  vertical 
distance  from  the  bottom. 

The  PSE  consists  of  a  bundle  of  thin-wall  PVC  tubing,  each  tube 
cut  to  a  precise  length  based  on  its  position  within  the  cross  section. 
The  length  of  each  tube  was  determined  in  the  following  way: 

Assuming  laminar  flow  through  all  the  tubes,  the  total  pressure 
drop  through  the  tubes,  P^,  is  given  by, 

Pj  =  Pf  +  Pe  A.2 

where  P^  is  the  pressure  loss  due  to  friction  and  P  is  the  pressure 

loss  due  to  entrance  and  exit  effects.  Px  and  P„  were  assumed  to  be 

T     e 

the  same  for  all  the  tubes,  while, 

P^  =  C  U(y,z)  L(y,z)  A. 3 
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where  C  is  a  constant  of  proportionality,  the  magnitude  of  which  depends 
on  the  fluid  viscosity  and  the  diameter  of  the  tubes.  Substituting 
equations  A. 3  and  A. 2  into  equation  A.l, 

at  z  =  D  and  y  =  0  equation  A. 4  becomes, 

^I^  =  l^max  LlO.f)  A.5 

Substituting  equation  A.5  back  into  A. 4: 

L(y,z)  =  L(0,D) . 

(2)1/7  (W::Z)V7  ^  g 

Time  and  labor  considerations  prevented  the  diameter  of  the  tubes 
from  being  very  small  even  though  the  requirements  of  laminar  flow 
dictated  the  opposite.  Thus,  tubes  with  1.91  cm  (0.75-inch)  I.D.  and 
2.13  cm  (0.84-inch)  O.D.  were  chosen.  The  tubes  were  stacked  in  a 
triangular  arrangement  alternating  between  28  and  27  tubes  per  row.  A 
total  of  32  rows  was  stacked  together  requiring  a  total  of  880  tubes. 

To  reduce  entrance  and  exit  effects  the  shortest  tube,  L(0,D),  was 
cut  nine  inches  long  so  that  all  the  tubes  would  have  a  length-to- 
diameter  ratio  greater  than  12.  The  length  of  each  tube  (in  inches) 
was  determined  by, 

L(y.z)  =     9.0 

(§)^/'  (^)^/'  A.7 
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The  requirement  that  the  flow  remain  laminar  determines  the  maximum 
allowable  velocity  through  the  PSE,  The  criteria  for  laminar  flow  are 
given  by, 

^  <  2200 

V 

2 

for  D  =  1.91  cm  and   =  0.01  ^ 

sec 

U    <  11 .5  cm/s 
max  — 

Due  to  the  many  assumptions  made  in  the  design  of  this  unit,  the 
resulting  velocity  profile  was  not  expected  to  match  the  equilibrium 
profile  exactly.  Indeed,  when  the  upper  layer  meets  the  salt  layer  a 
readjustment  in  the  velocity  profile  must  take  place.  The  intent  of 
the  PSE  was  to  minimize  longitudinal  variations;  it  was  not  expected  to 
eliminate  them  completely. 


APPENDIX  B 
INDIVIDUAL  DENSITY  PROFILES 

This  appendix  contains  the  individual  quasi-instantaneous  profiles 
of  density,  temperature,  and  Vaisala  frequency  obtained  by  the  constant 
speed  traverse  technique.  Throughout  the  appendix  there  are  two  plates 
for  each  run.  The  first  contains  the  profiles  of  density  and  tempera- 
ture, while  the  second  shows  the  corresponding  profiles  of  filtered 
density  and  Vaisala  frequency. 

The  vertical  coordinate  is  always  referenced  to  the  interface, 
defined  as  the  point  where  the  density  equals  the  average  density  of 
the  upper  and  lower  layers.  The  interface  is,  in  turn,  referenced  to 
the  bottom  of  the  tank.  The  density  is  given  in  terms  of  percent  spe- 
cific gravity  (S.G.  %),  defined  by: 

S.G.  %  =   (S.G.  -1.0)  100.0 
The  temperature  profile  is  plotted  as  the  deviation  from  the  mean 
temperature  (DELTA  TEMP.),  where  the  mean  temperature  is  the  average 
over  the  entire  profile.  The  filtered  specific  gravity  is  the  density 
after  the  least-square  smoothing  has  been  applied  (see  Chapter  III). 
The  Vaisala  frequency  is  nondimensionalized  by  the  maximum  value  of 
the  Vaisala  frequency  for  the  individual  profile,  NMAX. 
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APPENDIX  C 
THE  ERROR  IN  ASSUMING  Dp/Dt  =  0 

The  complete  continuity  equation  is  given  by. 


^  +  pv.tr  =0,  C.I 

where  the  material  derivative  is, 

D    9  ^  ^  ^ 
Dt  =  9t^  "•^' 

p  is  the  fluid  density  and  u  is  the  velocity.  The  density  of  salt 
water  is,  in  the  most  general  case,  dependent  on  the  specific  inter- 
nal energy,  e,  the  salinity,  S,  and  the  pressure,  P,  such  that 

D^  =  De  /9p_x      DS  / 3£>      DP  (dp_^  ^  ^ 

Dt   Dt  ^de'     ^   Dt  ^9S^^  .   Dt  ^3P^^  ^  ^'^ 

S,p        e,p        S,e 

The  equation  for  the  conservation  of  energy,  under  the  conditions 
that  there  is  no  internal  heat  source,  no  radiative  heat  transfer  and 
a  constant  coefficient  of  heat  conduction,  k,  is  given  by, 

De^^^.^^k^2^,y  C3 

Dt    p      p      p 

where  $  is  the  viscous  dissipation  function,  T  the  temperature,  and  y 
is  the  fluid  viscosity. 
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The  equation  governing  the  conservation  of  salt,  assuming  a  con- 
stant coefficient  of  diffusion,  is  given  by. 


DS    „2c  c^  -y 

P^  =  YgV  S  -  S^-u,  C.4 

where  j^   is  the  molecular  coefficient  of  diffusion  for  salt  in  water 
and  S  is  the  salinity  defined  as  the  mass  of  salt/unit  volume. 

Substituting  equations  C.2,  C.3,  C.4  into  equation  C.l  gives. 


^^  S,p  p^     ^^  S,p  p^ 

^SS-*^  n  ^^       3P  o  Q  f  Dt  ^-5 

e,P  e,S 


where    f  =   |^    ^+  (^) 

S,P  "^      e,p 


For  salt  water  at  atmospheric  pressure  and  room  temperature  the 

following  hold: 

3 

p  =  1 .0  g/cm 

k  ==  5.6  X  10  dynes/sec  °C 

-5   2 
Yg  =  1.1  X  10   cm  /sec 

P  -  10^  dynes/cm^ 

S  =  0.03  g/cm^ 

as   '-^ 

11=  4.3  X  10-''2g  sec^/cm^ 

If-  =  5.0  X  10"''^  sec^/cm^ 

AT  =  maximum  temperature  deviation  <  2°C 
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y$  is  the  energy  dissipation/unit  volume;  for  flow  in  a  channel 
the  following  estimation  is  made, 

% 

where  x  is  the  stress  on  the  walls  of  the  channel,  U   is  the  mean 
velocity  through  the  channel  and  R^  is  the  hydraulic  radius  of  the 
channel . 

For  the  flow  studied  herein, 

2 

T  <_  .1  dynes/cm 

U  £  10  cm/sec 

R,  =  40  cm 
h 

y$  £  ,03 

The  Kolmogorov  length  scale  which  corresponds  to  the  above  dissipation 
rate  is. 


X  =  ^       '.,.     =     .07  cm 
(.03)'/^ 

From  the  above  values  the  following  estimates  are  made: 


v^T  <  4  <  2^  ._     408 
A^  ~  (.07)^ 


v^S  <  4  <  "-^^  =  6.12 


^  <  pgU  <  10^ 
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f  =  4.3  X  10'^  +   .03  -  1.0  =  1.0 


and 


V-u  =   (4.3)(.03)(10"^^)  +  (4.3)(5.6)(4.08)(10"^) 
+   (1.0)(1.1)(6.12)(10"^)  +   (5.0)(10"^), 


therefore. 


->  ->  -4 

V-u  <  2  X  10     /sec. 


APPENDIX  D 
ENERGY  SPECTRA 

This  appendix  contains  the  energy  spectral  plots  of  the  horizon- 
tal and  vertical  velocity  fluctuations.  The  measurements  were  made 
during  four  of  the  nine  experiments  at  several  positions  above  (one  be- 
low) the  density  interface.  Two  plots  were  shown  for  each  run  nummber, 
the  first  is  the  energy  density  vs.  frequency  for  the  horizontal  veloc- 
ity flucturations  while  the  second  is  for  the  vertical  velocity  component. 

The  spectra  were  computed  by  a  fast  fourier  transform  (FFT)  algo- 
rithm using  the  subroutines  RHARM  and  HARM  from  the  Scientific  Sub- 
routine Package.  Before  the  spectra  were  computed,  the  mean  and  linear 

trend  were  removed  from  the  raw  data  by  linear  regression  and  the  ends 

2 
of  the  data  sets  were  tapered  using  a  sin  window  over  the  first  and 

last  10  percent  of  the  data  set.  The  data  sets  were  overlapped  50  per- 
cent in  the  time  domain.  Each  spectra  plot  shown  is  an  average  of 
several  individual  spectrum.  The  following  table  lists  the  pertinent 
facts  about  the  spectra: 


SAMPLE 

TNP^ 
16,384 

NSA^ 

EXP.    NO. 

DATE 

RATE 
20  Hz 

9 

12/6/78 

62 

3 

10/19/77 

30  Hz 

24,576 

93 

2 

5/3/77 

30  Hz 

24,576 

93 

1 

7/6/76 

20  Hz 

4,096 

15 

2Total  number  of  points  used  to  compute  spectrum. 
Number  of  spectra  averaged  to  obtain  spectral  plot. 
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APPENDIX  E 
VERTICAL  SHEAR  STRESS  PROFILE 

The  analysis  presented  in  this  appendix  is  carried  out  in  two  major 
steps.  First,  equations  2.1  through  2.4  are  developed  in  a  systematic 
way  from  the  differential  form  of  the  continuity  and  Navier-Stokes 
equations,  by  expanding  the  latter  equations  using  the  Reynolds  decompo- 
sition technique.  After  the  complete  equations  are  derived,  various 
terms  are  eliminated  until  equations  2.1  through  2.4  emerge.  This 
should  allow  the  reader  a  better  understanding  of  what  assumptions  are 
implied  by  the  use  of  equations  2.1  through  2.4. 

Following  this,  equations  2.1,  2.2  and  2.3  are  integrated  to 

determine  the  vertical  shear  stress  profile,  t  (x,y,z=0).  An  assump- 

xy 

tion  is  made  a  priori  to  simplify  the  integration.  The  velocity  and 
density  profiles  are  assumed  to  be  self-similar,  requiring  only  one 
length  scale,  6  and  6  respectively,  to  nondimensionalize  each  of  the 
profiles. 

The  continuity  equation  for  an  incompressible,  nondiffusive, 
nonconductive,  nondissipative  flow  is  given  by 

8x   ay   9z   ^'  ^' 

where  u,  v,  w  are  the  instantaneous  velocities  in  the  x,  y,  z  direc- 
tions, respectively.  In  a  real  fluid  equation  E.l  is  only  approximate- 
ly correct;  however,  for  the  range  of  conditions  studied  herein  the 
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error  resulting  from  the  use  of  equation  E.l  is  less  than  2  x  10"  /sec. 
Detailed  arguments  supporting  this  statement  are  presented  in 
Appendix  C. 

The  momentum  equations  for  an  incompressible,  Newtonian  fluid 
with  constant  viscosity  are: 


fau  ^  u  au  ^  V  au 
[dt  9x    ay 


W  3 


8V  ^  U  8V  ,  V  8V  ,  W  8V 

O  I  — r  +    +    +    

^  '  3t     3X     9y     3Z 


[3W  _!_  U 
3t 


9W  ^  V  3W  ^  W  3W 
3X     3y     3Z 


-3P 
3Z 


32w 


3X' 


32w 


3y^ 


3'^W 


32^ 


E.2 


E.3 


E.4 


where  p  and  P  are  the  instantaneous  density  and  pressure,  respectively, 
y  is  the  fluid  viscosity  and  g  is  the  gravitational  acceleration.  The 
variation  in  the  fluid  viscosity  resulting  from  salinity  and  tempera- 
ture variations  is  less  than  four  percent  for  the  flow  studied  herein, 
therefore,  the  fluid  viscosity  is  considered  constant. 

Employing  Reynolds  decomposition  technique,  the  dependent 
variables  are  represented  as  the  sum  of  fluctuating  and  time 
independent  components  as  follows: 


p(x,y,z,t)  =  po(x,y,z)  +  p'(x,y,z,t), 

P(x,y,z,t)  =  Po{x,y,z)  +  P'(x,y,z,t), 

u(x,y,z,t)  =  U(x,y,z)  +  u'(x,y,z,t), 

v(x,y,z,t)  =  V(x,y,z)  +  v'(x,y,z,t), 

w(x,y,z,t)  =  W{x,y,z)  +  w'(x,y,z,t). 


"N 


> 


E.5 


J 
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The  time  average  of  any  fluctuating  quantity,  f'(t),  is  defined  by. 


«+T 

f'(t)  =  l™'lzr  J  _^  ^'(*)^*' 


and  is  denoted  in  the  equations  to  follow  by  an  overbar.  By  definition. 


p'  =  P'  =  u'  =  v'  =  w'  =  0. 
Substituting  equations  E.5  into  equation  E.l  gives 


3U^3V^9W^8u'^3v'^3w'^f,  p, 

3x   ay   8z   3x    dy         az  * 


Time  averaging  equation  E.6  reduces  it  to 


ax   ay   az    '  '^•' 


and  subtracting  equation  E.7  from  equation  E.6  results  in 

3u^^  ^  ^  awl 

ax   ay   az  "^'^ 

Thus,  not  only  is  the  dilatation  of  the  mean  flow  zero,  but  at  any 
instant  in  time,  the  divergence  of  the  turbulent  velocity  fluctuations 
is  also  zero. 

In  a  parallel  manner,  the  momentum  equations  governing  the  mean 
flow  of  a  turbulent  stratified  fluid  are  derived  by  substituting 
equation  E.5  into  equations  E.2,  E.3  and  E.4.  After  time  averaging 
and  much  simple  but  tedious  manipulation  using  equations  E.7  and  E.8, 
equations  E.2,  E.3  and  E,4  become: 
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1   + 


p'u 


PoU 


9X 


1    + 


P'V 

PoV 


3y 


1  + 


p'w' 

PoW 


«l^ 


^U_3_p'u'  +V_a_p'u'  +W_3_p'u' 

Po    3X  po    9y  Po    9Z 


^l_^p'u'u'   ^]_^p'v'u'   ^l_3_p'w'u' 
Po    9X  Po    9y  Po    9Z 


=  zilPo.  +    3_  -l^'lJ'    +  1!_  3^U    _    U'u'     9po 
Po    9X  9X  po    9X^    ~      Po      9X 


_^    9_  -U'V'    ^   y_  9^1] 

3y  Po  "9y^ 


u'v'    9po 
Po     9y 


^    9_  -U'w'    ^   M      9^U 


U'W'     9po 


9Z 


Po    9Z^  Po      9Z 


E.9 


1  + 


p'u' 

PoU 


uf  + 

9X 


1    + 


p'V 

PoV 


ay 


1  + 


p'w' 

PoW 


w|^ 


+  U_9_p'v'   +V_9_p'v'   +W_i_p'v' 
Po    9X  Po    9y  Po    3Z 


^l_l_p'u'v'   +l_9_p'v'v'   +T_8_p'w'v' 
Po    9X  Po    9y  Po    9Z 


•1    9Po    _+^-u'v'^y_  9£V   _    u'v'    9p 


Po     9y 


9X 


Po    9X  Po      9X 


^9_-v'v'    ^  1L_  9^V         V'v'     9pc 

9y  po  9y^  "    po    9y 


_,    9      -W'v'    ^   M      92V        W'V'    9po 


9Z 


Po    9Z^  Po      9Z 


E.IO 
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1  + 


PoU 


1  + 


P'V 

PoV 


3y 


1  + 


p'W 

PoW 


W  9W 
"  3Z 


^U_3_p'w'  +  V_i_p'w'  +W_l_p'w' 

po  9X        po  9y        Po  3Z 


^  l_a_p'u'w'  ^  l_3_p'v'w'  ^  l_3_p'w'w' 

Po  9X  Po  9y  Po  32 


=  ll  9Po.  +  1_  -U'W  ^  ;j_  9^W  _  U'w'  dp_ 


Po  3Z 


9X 


Po  9X 


7" 


Po   3X 


.  9   -V'w'  ,  y   g^W   V'W'  9po 

+  - —         + : — 17- ■ 


9y 


Po  9y^ 


Po  3y 


_^  9_  -w'w'  _^  jj_  9^W  -w'w'  dpj 
3Z         Po  3Z^       3Z 


E.11 


Equations  E.7  through  E.ll  have  been  developed  from  the  differen- 
tial form  of  the  continuity  and  momentum  equations  while  invoking  a 
minimum  number  of  assumptions;  therefore,  these  equations  enjoy  gen- 
eral applicability.  However,  because  of  their  complexity  and  the 
problem  of  closure,  the  system  cannot  be  solved.  In  order  to  reduce 
the  equations  to  a  form  which  can  be  solved,  assumptions  must  be  made 
regarding  the  nature  of  the  specific  flow  of  interest. 

First,  if  at  the  center  of  the  tank  (z=0),  where  the  measurements 
presented  in  Chapter  IV  were  made,  the  flow  is  assumed  to  be  nearly  two 
dimensional  and  free  of  any  secondary  currents,  then 


U  _  9W  _  9U  _  3V  _  9Po  _  9po  _  -T-TT  _ 


9Z    3Z    9Z 


9Z 


3Z 


0 
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Under  these  assumptions,  equation  E.7  reduces  to 

Next,  if  the  pressure  throughout  the  fluid  is  assumed  to  be  hydro- 
static, equation  E.IO  reduces  to 

ly^  =  -g  Po(x,y)  E.13 

Finally,  if  the  mean  advective  acceleration,  the  horizontal  pres- 
sure gradient,  the  side  wall  shear  stress  and  the  interfacial  shear 
stress  are  assumed  to  be  the  only  significant  terms  in  the  x-component 
momentum  equation  and  if  po  is  approximately  equal  to  p  ,  so  that  the 
Boussinesq  approximation  holds,  then  equation  E.9  becomes 

^  +  lUV  ^  ^  3Po  ^  9_  T     8_  T 

3x    9y   p^   ax    9y  ^^       9z  ^^'  "-'^ 

where 


XZ  p^  3Z 

is  the  shear  stress  produced  by  the  side  walls  and 


xy        p^^  3y 

is  the  shear  stress  in  the  plane  z=0. 

Equation  E.14  can  be  further  simplified  if  the  flow  does  not  change 
with  the  distance  downstream  and  the  side  wall  stresses  are  insignificant. 
Under  these  conditions  equation  E.14  becomes 

-1^=P^.  E.15 

p  9x   ay 
u      -^ 

For  an  ideal  two  layer  freesurface  flow,  equations  E.13  and  E.15 

predict  a  linear  shear  stress  profile  which  is  zero  at  the  freesurface 

and  maximum  at  the  interface.  The  more  complete  model,  given  by 
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equations  E.12  through  E.14,  accounts  for  the  side  wall  friction, 
nonhomogenous  density  and  advective  accelerations  due  to  boundary  layer 
development.  Next,  these  equations  are  solved  to  determine  the  shear 
stress  profile  at  the  center  of  the  tank  (z=0). 

Intergrating  equation  E.12  from  the  free  surface,  y  =  v(x),  to  y 
gives 


v(x) 
V(x.y,0)  =  /  1^  U(x,y)dy  +  U(x,v)  g  ,  E.16 


/ 


y 

under  the  assumption  that  a  fluid  particle  at  the  surface  remains  at  the 
surface  (i.e.  V{x,v)  =  Ll(x,v)  dv/dx). 

Invoking  the  self-similarity  hypothesis,  the  mean  velocity  profile 
can  be  defined  by 

U(x,C*)  ^  U(x)  f(c*),  E.17 

where  f  is  a  nondimensional  function  of  the  nondimensional  variable 
C*,  U(x)  is  the  freestream  velocity,  and 

'*  =  ^-n(x)  ^  y-D(x)-n(x) 


^    V^)     V"^ '  E.18 


where  D(x)  is  the  distance  from  the  bottom  of  the  tank  to  the  density 
interface,  n(x)  is  the  vertical  shift  between  the  density  and  velocity 
profiles  and  6^(x)  is  the  thickness  of  the  velocity  transition  layer 
above  the  interface.  By  definition,  in  the  freestream, 

fU*)  =   1   ;  ^^>  C*  >  2  . 

u 

Making  the  proper  substitutions,  equation  E.16  can  be  rewritten  as 

(v-D-n)/6^ 

V{x,c*,0)  =  6^(x)  /  |-u(x,?*)d5*  +  U(x)  ^  ,  E.19 


203 


where 


and 


1^  U(x,5*)  =  ^^  f  (?*)  +  U(x)  f  (5*)  If  ,  E.20 


9|j^  _  -5*   ^   JL  ^  dD  .  dn  X  CO. 

dx        J^)   dx  "6  ^  dx   dx  ^  •  ''•'^' 


Substituting  equations  E.20  and  E.21  into  E,19  gives 


■/ 


(v-D-n)/6^j  (v-D-n)/6^ 


V(x,5*,0)  =  6^(x)  ^)  /f(c*)  dc*  -  U(x)  ^/  f  (S*)  5*  d?* 


-U(x)  (f +^)  (  1  -f(C*)  )+U(x)^. 

Integrating  the  second  term  on  the  right  hand  side  by  parts  gives 

/(v-D-n)/6^j 
f(?*)  d5* 


u^"'  dx    ^^'^'   dx 

5* 


+  U(x)  ^  5*f(C*)  -  U(x)  (  f  +  ^  )(  1  -  f(C*)  ) 
^  ^(x)  ^  -  U(x)  ^  (^^)  .. 

Integrating  equation  E.13  from  the  freesurface,  y  =  v(x),  to  y  gives 

v(x) 
Po(x,y)  =  g  /  Po(x,y)dy  +  Poi^v)  .  E.23 


■/ 
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Differentiating  equation  E.23  with  respect  to  x,  gives,  after  applying 
Leibnitz's  Rule, 

v(x) 


•!'■ 


|^>  =  g  t^y  ^  9  P»(x,v)  f  . 


y 

If  the  density  distribution  is  also  self-similar,  it  can  be  represented 
by 

Po(c"^)  =  p^j  +  Ap  Q*U'^),  E.25 

where  p*  is  a  nondimensional  function  of  the  nondimensional  vertical 
coordinate  t,  ,   Ap  is  the  density  change  across  the  interface  and 

r"^  -  f   -  y-D(x)  F  oc 

where  D(x)  is  the  distance  from  the  bottom  of  the  tank  to  the  density 
interface  and  6(x)  is  the  thickness  of  the  density  interface.  By  defi- 
nition 

p*  =  0   ;  5"^  >  3 

P*  =  1   ;  c""  1  -3 
Making  the  above  substitutions,  equation  E.24  becomes 

(v-D)/6 


9Po(x,g"^)  _  ^^   d6  ,  ^ 

+ 


^P*U^)   +    P*U^)  dg^  ) 


^9AP^P*U^)^gp„S.  E.27 

Next,  assuming  the  shear  stress  maintained  by  the  side  walls  can  be 
written  as  a  linear  function  of  z  times  the  value  of  the  stress  at  the 
side  wall ,  then 
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^xz^^'-^'^^  =  -  T^(x,y)  z/W,  E.28 

where  W  is  the  half  width  of  the  tank  and  x  is  the  shear  stress  at  the 

w 

side  walls.  Differentiating  equation  E.28  with  respect  to  z  gives, 

Nz  _   S/^) 

9Z     "    W 

Furthermore,  if  the  side  wall  shear  stress  can  be  represented  by 
Tjx,y)  =  Mx^)  U2(x,y,0)  =  ^-%^)  ^(x)  f^^*)    , 


where  fo  is  the  side  wall  friction  factor,  then 


^=  "  ^iw'^*^  U^(x)  ^^U*)-  E.29 


The  value  of  fo  depends  on  the  Reynolds  number  and  may  depend  on  the 
density  distribution.  Therefore,  fo  may  be  a  function  of  5*. 

Integrating  equation  E.14  from  the  free  surface  to  y  gives, 

v(x) 

1^  U^(x,y)dy  +  U(x,v)  V(x,v)  -  U(x,y)  V(x,y) 


^  '  ^^°   dy  +  T„„/,v)  -  T, Jx,y)  +  I  I-  T,  fx,y)  dy. 


p   /  9x   -^    xy  '  '    xy^    '-^ '       J     9x  xz 

y  y 

The  above  equation  can  be  transformed  into  the  c*  and  5  coordinate 

system  by  substituting  equations  E.17,  E.18,  E.21,  E.22,  E.27,  and  E.29. 
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The  resulting  equation  is 


T^^(x,c*,5'')  =  -  g  ^  ^(^  -  C"")    (Slope  of  freesurface)    E.30 

(v-D)/6 

-  g^^  6/  p*(c'^)dc'^    (Slope  of  interface) 
X  j^ 

(v-D)/6  (V-D)/6 

(Diffusion  of  I.F.)  -  g  f-^  5[   /  pMs"")  dcV  +  h"'p*(c"')  d?""  ] 


(v-D-n)6^ 
(Side  wall  shear)   -  5^)  6„   /  fof^(5*)  dc* 


■J' 


(Advective  acceleration  +  LJ^(x)  (  ^  +  ^  "  ^  )  (  1  "  ^(s*)  ) 
due  to  nonuniformity  (v-D-n)i5        (v-D-ri)(5 

of  the  flow)  +U^&^lfU*)lfU*)6g*-2      /f^(?*)dc*] 


)  jfU*)   6g*  -  2      If 


(Advective  acceleration  Jv-D-n)/6     (v-D-n)/5 

due  to  growth  of  the    +  U'^(x)  j-^  [f (c*)  /  f  (5*)  dc*  -     f  U*)   d?* 


boundary  layer) 


+  (^)  (1  -  f(c*)  )]. 
u 


Given  the  velocity  and  density  profile  at  a  fixed  value  of  X, 
equation  E.30  becomes  a  linear  equation  in  seven  parameters. 


207 


The  seven  parameters  are: 

•g-   the  slope  of  the  freesurface, 

•^   the  slope  of  the  interface, 

^   the  horizontal  gradient  of  the  mean  velocity  in  the  uniform 
portion  of  the  velocity  profile. 


the  spacial  rate  of  growth  of  the  velocity  boundary  layer. 


d5 

u 

dx 

fo   the  side  wall  friction  factor. 


■^      the  spacial  rate  of  growth  of  the  separation  distance  between  the 
density  and  velocity  interfaces, 

-^      the  spacial  rate  of  growth  of  the  density  transition  region  due 
to  diffusion. 

The  values  of  the  parameters  are  determined  as  follows: 

fo   was  set  equal  to  .02  (the  value  for  hydraulical ly  smooth  walls). 
Several  values  ranging  from  .01  to  .035  were  tried  and  the  effect 
on  the  shear  stress  profile  was  always  minor. 

-rr      was  determined  from  experimental  data  to  be  about  2.5  x  10"  (see 
the  section  in  this  chapter  on  the  density  structure).  Again, 
several  values  ranging  from  0  to  10-3  were  tried  with  little 
effect  on  the  shear  stress  profile. 

■j-  and  J—   The  data,  as  well  as  visual  observation,  indicates  that 
the  velocity  drops  to  a  yery   small  value  or  zero  at  a 
point  1  to  3  centimeters  below  the  interface.  Since  there 
is  no  motion  below  this  point,  the  equations  of  motion  re- 
duce to  a  balance  between  the  hydrostatic  pressure  gradients 
imposed  by  the  freesurface  set-down  and  the  interfacial 
set-up.  dv/dx  and  dD/dx  were  determined  by  requiring  the 
shear  stress  to  be  zero  in  the  region  of  no  motion. 

dU  d   ^"^u 

dx'  dx'  dx~   These  parameters  were  determined  by  either  fitting  the 

shear  stress  to  the  experimental  data  using  a  least 

square  error  technique  or  by  fixing  one  or  more  of  the 

three  parameters  and  solving  for  the  remaining  parameters 

by  the  least  square  approach. 

The  results  of  this  analysis  are  shown  in  Figure  33. 
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